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Abstract. A methodology for precise determination 
of the fundamental geodetic parameter w 0 , the po¬ 
tential value of Gauss-Listing geoid, as well as its 
time derivative w 0 is presented. The method is 
based on: (i) ellipsoidal harmonic expansion of the 
external gravitational field of the Earth to de¬ 
gree/order 360/360 (130321 coefficients 

http://www.uni-stuttgart.de/gi/research/index.html- 
#Projects) with respect to the International Refer¬ 
ence Ellipsoid WGD2000, at the GPS positioned 
stations, (ii) ellipsoidal free-air gravity reduction of 
degree/order 360/360 , based on orthometric 
heights of the GPS positioned stations. The method 
has been numerically tested for the data of three 
GPS campaigns of Baltic Sea Level Projects (epoch 
1990.8, 1993.4, and 1997.4). New w 0 and w 0 val¬ 
ues (w 0 = (62636855.75 ± 0.21) (m 2 / s 2 ) , w 0 = 
(-0.0099 ± 0.00079) (m 2 / s 2 ) / year , u ; 0 !l = 

6379781.502(ra), w 0 / 7 = 1.0(mm / year ), and 
7 = —9.81802523m / s 2 ) for the test region (Bal¬ 
tic Sea) have been obtained. As by-products of the 
main study we also succeeded (i) to determine the 
high-resolution Sea Surface Topography map for 
the Baltic Sea, (ii) to determine the most accurate 
regional geoid amongst four different regional 
Gauss-Listing geoids currently proposed for the 
Baltic Sea , and (iii) to determine the differences 
between National Height Datum of countries 
around the Baltic Sea. 

Keywords: Potential value of geoid. Time varia¬ 
tions of the potential value of geoid. Sea Surface 
Topography. Differences between National Height 
Datums. 


0 Introduction 

There are four quantities, namely {w 0 ,gm,j 2 ,cu}, 

which are currently accepted by geodetic commu¬ 
nity, as the fundamental geodetic parameters, as 
documented by E. Groten (2000), for instance. 


Among these parameters, w 0 plays a crucial role in 
the geoid determination as well as the computation 
of the best fitting reference equipotential surface to 
the geoid. For example, the reference potential 
field W(r) — gm / r = w 0 gauged to the geoid 

potential value w 0 , leads to the sphere S% =gm / WQ as 

the best fitting reference equipotential surface to 
geoid , or the Somigliana-Pizzetti potential field 
gauged to the geoid potential value w 0 generates 

the ellipsoid of revolution £^ a b , i.e. WGD2000 (E. 

Grafarend and A. Ardalan 1999), as the best fitting 
reference equipotential surface to the geoid. Ac¬ 
cording to E. Grafarend and A. Ardalan (1999, 
page 614), the accuracy of such a best-fitting ellip¬ 
soid is mainly driven by the accuracy of w 0 . 

Here we are proposing the following algorithm 
for the precise determination of fundamental geo¬ 
detic parameter w 0 : 

(i) determine the gravity potential at the accurately 
positioned GPS stations, by using the 
ellipsoidal harmonic expansion of degree/order 
360/360 (130321 coefficients http://www.uni- 
stuttgart.de/gi/research/index.html#Projects), 
plus the ellipsoidal centrifugal potential. 

(ii) reduce the derived potential values to the Mean 
Sea Level (MSL) via the precise ortho-metric 
height of GPS positioned stations, and the 
ellipsoidal free air gravity reduction of de¬ 
gree/order 360/360. 

The input data for the “ w 0 operational procedure” 
quoted above can be obtained by GPS observations 
at permanent tide gauges or at stations which are 
located close to tide gauge stations and are well 
connected to the tide gauge stations by precise lev¬ 
elling. The GPS observations collected for Baltic 
Sea Level Project (see for instance, J. Kakkuri, 
(1990, 1995) and M. Poutanen and J. Kakkuri , 
(1999)) provide us with such a source of informa¬ 
tion. Figure 0-1 presents a typical configuration of 
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the GPS and permanent tide gauge stations of Baltic 
Sea Level Project. 

Intuitively speaking, the global geopotential 
model together with position information derived 
by GPS provide us with the global gravity informa¬ 
tion, and, in contrast, the connection of GPS sta¬ 
tions to the tide gauge by precise levelling, supplies 
us with the Mean Sea Level information. Recalling 
that the geoid according to the Gauss / Listing defi¬ 
nition (C. Gauss 1828 page 49, J. Listing 1873 page 
45) is an equipotential surface which fits to the 
Mean Sea Level in an optimal way, the approach 


outlined above should be adequate for a global 
evaluation of w 0 . However, since the tide gauge 
observations, aside from random errors, are always 
affected by some local effects, we may obtain dif¬ 
ferent w 0 values for different tide gauge stations. 
However, the average value of w 0i data (the index 
i runs from one to the total number of the tide 
gauge stations) computed from a set of globally 
well distributed tide gauges may be considered an 
unbiased or, at least, a minimum biased estimate of 
the geoid potential value w 0 


GPS station 


measured height 
via precise levelling 



Figure 0-1: The GPS station, permanent tide gauge station, and measured precise or¬ 
thometric height of the GPS station above Mean Sea Level. 


In Section 7, we set up the theoretical back¬ 
ground of our w 0 computation methodology. In 
Section 2 we shall introduce those data that we will 
use for various numerical computations. Sections 3 , 
4 and 5 focus on the results of w 0 computation 
based on GPS data of Baltic Sea Level Project , 1 st 
campaign, 2 nd campaign, and 3 rd campaign. Section 
5, leads to the best estimate of w Q = 

(62636855.75 ± 0.21)(ra 2 /s 2 ) with a remarka¬ 
bly small mean square error. In Section 6 , we aim at 
an estimate of the time derivative of w 0 , namely 

w Q = (0.0099 ± 0.00079) (m^/s^)/year. Section 7 
discusses the apparent difference in the height da¬ 
tum between various countries around the Baltic 
Sea as well as the transfer function of gravity po¬ 
tential differences into geometric heights in order to 
obtain the differences between height datums. For 
the unification of national, regional or continental 
to a global height datum we acknowledge contribu¬ 
tions of R. Rummel and P. Teunissen (1988), B. 
Heck, and R. Rummel (1990), P. Xu and R. Rummel 


(1991), P. Xu (1992), R. Rapp (1994), R. Rummel 
and K. Ilk (1995), and F. Sansd and S. Usai (1995). 
Finally Section 8 concentrates on the determination 
of Sea Surface Topography of Baltic Sea at the tide 
gauges of the Baltic Sea Level Project. 

1. W 0 / W 0i Computations: 

Problem Formulation 

It is well known that a slowly uniformly rotating, 
self-gravitating liquid body of radial mass distribu¬ 
tion, forms in its equilibrium state an ellipsoid of 
revolution at its boundary. Therefore, the Earth as a 
massive body, which once was a liquid at hydro¬ 
static equilibrium, globally does not deviate consid¬ 
erably from an oblate spheroid or ellipsoid of 
revolution. Therefore, the external gravity field of 
the Earth can most conveniently be presented in 
terms of an ellipsoidal model as explained in Table 
1 - 1 . 
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Table 1-1: Additive decomposition of the gravity potential W into gravitational potential U and the cen- 
trifugal potential V . _ 


where 


and 
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According to Table 7-7, we have additively de¬ 
composed the gravity potential W into the gravita¬ 
tional potential U and the centrifugal potential V ; 
both represented in terms of Jacobi ellipsoidal co¬ 
ordinates {A, cf),u} also called spheroidal. In the 

space external to the Earth, its gravitational poten¬ 
tial is harmonic , namely it is an element of the 
three-dimensional Laplace equation with associ¬ 
ated Legendre functions of the first Pnm (sin <j>) and 

of the second kind Q* nm ( u / b ) as well as { cos mA, 
sin | m | A } as eigenfunctions. Accordingly, in 
terms of surface spheroidal harmonics e n m(A,^>) 
given by (1.4), the equations (1.1)—(1.3) are proper 
representations. For the definition of the Jacobi el¬ 
lipsoidal coordinates {A, </>, u} we refer to Appen¬ 
dix A. Appendix B provides us with the definition 
of the normalised associated Legendre functions of 
the first as well as of the second kind. Finally, Ap¬ 
pendix C outlines how the ellipsoidal harmonic co¬ 


efficients Unm can be obtained from the available 
spherical harmonic coefficients. 

Equation (1.1) is exact and without any ap¬ 
proximation. However, in practice we have to ap¬ 
proximate it by using a limited number of ellipsoi¬ 
dal harmonic coefficients, say up to degree/order 
360/360. Consequently, we arrive at an approxi¬ 
mate form of (1.1), which can serve as a reference 
gravity field of the Earth. Such a reference field 
covers the gravity effects of the features which are 
bigger than 50-60km. 

Roughly speaking, now if we place the GPS re¬ 
ceiver at point p Q on the Mean Sea Level (i.e. 
geoid) and measure the Cartesian coordinates 
{x 0 ,y Q ,z Q } of p 0 , then once we have converted 
these Cartesian coordinates into Jacobi ellipsoidal 
coordinates {Aq,<A),^ 0 } (via the transformation re¬ 
lations given in Appendix A) an estimation of geoid 
potential value w 0 can be obtained according to 
Table 1-2. 


Table 1-2: High degree/order reference potential field of the external gravity field of the Earth. 
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Since at the geoid’s level the geometrical varia¬ 
tion of the equipotential is quite smooth, de¬ 
gree/order 360/360 is enough for the computation 
of the geoid potential value at the level of accuracy 


better than 0.5m 2 / s 2 . Note that the high- 

frequency changes of the gravity field is mainly due 
to the topographical masses, therefore at the geoid’s 
level we do not have any high-frequency variations 
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and as such the ellipsoidal harmonic expansion of 
degree/order 360/360 can satisfactorily represent 
the geoid potential value w 0 . 

In practice however, we may place the GPS re¬ 
ceive at point p , which is at the vicinity of a tide 
gauge station, and measure the Cartesian coordi¬ 
nates {x,y,z}. These Cartesian coordinates once 

converted into Jacobi ellipsoidal coordinates 
{A, </>, u} can provide us with the gravity potential 


W{\(j),u) at the point p, which is of course dif¬ 
ferent from the geoid potential value w 0 . 

Knowing, that the potential difference between 
the points p{\<j),u} and the point Po{\?<A>>^o} 
on the geoid is only due to the orthometric height of 
p{ A, </>, u} , the geoid potential value w 0 can be 
written in terms of the Taylor series expansion 
around the potential value W(\(f),u) at the point 

p{ A, </>, u} as explained in Table 1-3 . 


Table 1-3: Taylor series expansion of the geoid potential value w 0 around the potential value W{\,(j),u) at 
the point p{ A, </>, u} . 

“Taylor series expansion of the geoid potential value w 0 ” 
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In Table 1-3 , D U W is the partial derivative of 
potential W(\,(f>,u) with respect to u, and 
(u 0 — u ) is the height of GPS station above MSL, 

of course in terms of Jacobi ellipsoidal height com¬ 
ponent. The partial derivative D U W is related to di¬ 
rectional derivative of W along the coordinate line 
of u as introduced by (1.7). For the definition of 
directional derivative V e „ W along the coordinate 
line of u, we refer to equation (1.8). 
e A = t a/V&T, = 1and 

e u = t u / are orthonormal base vectors of the 
Jacobi ellipsoidal coordinates {A, </>, u }, and Jiff » 
^fg~ , and are metric tensors. See Appendix 
A, equation (A. 10) for the definition of the metric 


tensors in terms of Jacobi ellipsoidal coordinates 
{A, </>, u} . The vectors t A = dx / d\ , 
= dx / dcf) , t u = dx/du are tangent to coor¬ 
dinate lines of Jacobi ellipsoidal coordinates A , <fi , 
and u, respectively. Therefore, (1.6) up to the 
terms of the order of magnitude O((u 0 — uf) can 
be written as (1.9). 

According to C. Eringen (1962, page 437) 
= V Quu (uq — u) is the physical component 
of the Jacobi ellipsoidal height difference (u 0 — u ). 

Partial derivative of associated Legendre functions 
of second kind with respect to u can be derived 
from the recursive relations of Table 1-4. 


Table 1-4: Partial derivative of associated Legendre functions of second kind with respect to u. 
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We call the second term of (1.9), i.e. 
A W (A, (j), u ), the “ellipsoidal free-air reduction ” 

of the potential value W(\(j>,u) to the geoid’s sur¬ 
face. Since in the ellipsoidal harmonic expansion all 
masses of the Earth are condensed inside the refer¬ 
ence ellipsoid of WGD2000, in the eyes of ellipsoi¬ 
dal harmonic model, we have actually free-air be¬ 
tween p and p 0 ! The physical component 


■fg^(u 0 — u) here can be interpreted as the or¬ 
thometric height H 0 of the GPS station with oppo¬ 
site sign, i.e. — u) = -H 0 . 

Therefore, the geoid potential value w 0 can be 
derived from the Jacobi coordinates {A, 0, u} of the 
GPS station p and the orthometric height of p, 
, as outlined in Table 1-5. 
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Table 1-5: Operational procedure for computing the geoid potential value w 0 
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Having set up the theoretical foundation of the w 0 
computation problem, we can begin our case study 
by the computing w 0 and w 0 from the GPS obser¬ 
vation of the Baltic Sea Level projects, first, sec¬ 
ond, and third campaigns in the next sections. 

2. W 0 Computations: Input Data 

Based on the method described in previous section 
here we compute the w 0 value via the GPS obser¬ 
vations of three successive GPS campaigns of Bal¬ 
tic Sea Level Project. For a review of the state of 
the art of Baltic Sea Level project J. Kakkuri (1990, 
1995), and M. Poutanen and J Kakkuri (1999) are 
recommended. Table 2-1-Table 2-3 present the 
Cartesian coordinates of the GPS stations of Baltic 
Sea level project, 1 st , 2 nd , and 3 rd campaigns, respec¬ 
tively. Table 2-4-Table 2-6 are presenting the 
computed Jacobi ellipsoidal coordinates (see Ap¬ 
pendix A for transformation relations) of those sta¬ 
tions given in Table 2-1-Table 2-3. The ellipsoidal 
harmonic coefficients needed for the series expan¬ 
sion of (1.15) and (1.16) are provided via the trans¬ 
formation of the spherical harmonic coefficients of 
EGM96 ( F. Lemoine et al 1998) to ellipsoidal ones 
with respect to the reference ellipsoid of WGD2000 


through exact transformation relations given in Ap¬ 
pendix C. The ellipsoidal harmonic coefficients are 
computed in mean-tide permanent-tide system (see 
M. Ekman (1996) for the definition of various per¬ 
manent-tide systems). For this purpose we first 
transferred the second zonal spherical harmonic co¬ 
efficients of EGM96 geopotential model from tide 
free system into mean tide system, and then applied 
the transformation machinery of Appendix C to ob¬ 
tain ellipsoidal harmonic coefficients from spherical 
harmonic coefficients. We have made available the 
computed ellipsoidal harmonic coefficients to pub¬ 
lic, and interested readers can download the ellip¬ 
soidal harmonic coefficients plus the manual of 
using them, from the homepage of Geodetic Insti¬ 
tute of the Stuttgart University (http: // www.uni- 
stuttgart.de / gi / research / index.html # Projects). 
As the reference ellipsoid we use the WGD2000 in 
mean tide system (c.f. E. Grafarend and A. Arda¬ 
lan , 1999). In other words, we performed all the 
computations in mean-tide permanent-tide system. 

We obtained the orthometric heights, of GPS 
stations from five different sources namely: M. 
Poutanen et al. (1999), J. Kakkuri (1995), J. Kak¬ 
kuri and M. Poutanen (1997), M. Poutanen et al. 
(1999), and J. Kakkuri (2000). Consequently, we 
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performed the computations under five cases de¬ 
fined below. 

Case 1: the orthometric heights from 
M. Poutanen et al. (1999) 

Case 2: the orthometric heights 
from J. Kakkuri (1995) 

Case 3: the orthometric heights 

from J. Kakkuri and M. Poutanen (1997) 
Case 4: the orthometric heights 

from M. Poutanen et al (1999) 

Case 5: the heights from Kakkuri (2000) 

In case 1, the orthometric heights H 0 are given in 
their respective national height systems, measured 
directly by precise levelling to tide gauges, while 
for the cases 2-5 orthogonal heights H 0 are de¬ 


rived from GPS ellipsoidal height h and different 
geoid solutions N proposed for Baltic Sea. Recall 
that neglecting the curvature of plumb line follow¬ 
ing relation holds. 

H 0 =h-N (2.1) 

It should be mention that the difference between 
orthometric heights of cases 2-5 is mainly due to 
different geoid solutions and application of differ¬ 
ent reductions to obtain mean geoid via polynomial 
fitting, and using different starting latitudes (see for 
example M. Vermmer (1995) for more details). 

In the following sections, the results of w 0 compu¬ 
tations for each of the above mentioned cases will 
be presented. 


Table 2-1: Cartesian coordinates of GPS stations of BSL project first campaign 1990.8, in ITRF 91 reference 
frame. 


Station 

Name 

X 

(m) 

Y 

(m) 

Z 

(m) 

Borkum (Ger) 

3770668.4100 

446076.5560 

5107686.4450 

Degerby (Fin) 

2994005.5910 

1112565.2290 

5502270.8220 

Furuogrund (Swe) 

2527022.8530 

981957.3370 

5753940.7500 

Hamina (Fin) 

2795471.3860 

1435427.7410 

5531682.1280 

Hanko (Fin) 

2959173.1020 

1254706.6550 

5490604.5210 

Helgoland (Ger) 

3706045.0350 

513713.1770 

5148193.1900 

Helsinki (Fin) 

2885134.8290 

1342693.6440 

5509043.9240 

Kemi (Fin) 

2397170.2880 

1093246.9200 

5789077.1960 

Klagshamn (Swe) 

3527585.8260 

807513.8050 

5234549.4560 

Klaipeda (Lit) 

- 


- 

Kronstadt (Rus) 

- 

- 

- 

List/Sylt (Ger) 

3625340.1590 

5378854.8410 

5202539.3040 

Mantyluoto (Fin) 

2831096.8750 

1113102.7630 

5587164.9240 

Molas (Lit) 

- 

- 

- 

Olands N. U. (Swe) 

3295551.6790 

1012564.8840 

5348113.5190 

Raahe (Fin) 

2492699.6820 

1131503.6280 

5741504.1800 

Ratan (Swe) 

2620087.7400 

1000008.3720 

5709322.5000 

Shepelevo (Rus) 

- 

- 

- 

Spikama (Swe) 

2828573.5470 

893623.7640 

5627446.8610 

Stockholm (Swe) 

3101011.4900 

1013009.1090 

5462375.0730 

Swinoujscie (Pol) 

3649458.4740 

927709.8850 

5130741.4690 

Ustka (Pol) 

3545014.4560 

1073939.6850 

5174949.7730 

Vaasa (Fin) 

2691307.2790 

1063691.5640 

5664806.2010 

Visby (Swe) 

3249304.5320 

1073624.8010 

5364362.8610 

Warnemiinde (Ger) 

3658230.7910 

783507.3140 

5148395.8730 

Ger: Germany, Fin: Finland, Swe: Sweden, Lit: Lithuania, Pol: Poland, Rus: Russia 
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Table 2—2: Cartesian coordinates of GPS stations of BSL project second campaign 1993.4, in ITRF 93 refer¬ 
ence frame. 

Station X Y Z 

Name _(m)_(m)_(m)_ 

Borkum (Ger) 3770668.1000 446076.4240 5107686.2360 

Degerby (Fin) 2994005.4960 1112565.2460 5502270.9240 

Furuogrund (Swe) 2527022.9170 981957.2520 5753940.9600 

Hamina (Fin) 2795471.2440 1435427.7570 5531682.2020 

Hanko (Fin) 2959172.9990 1254706.6680 5490604.6630 

Helgoland (Ger) 3706044.9580 513713.1520 5148193.3740 

Helsinki (Fin) 2885134.7810 1342693.6430 5509044.0630 

Kemi (Fin) 2397170.2290 1093246.8740 5789077.1460 

Klagshamn (Swe) 3527585.8220 807513.8260 5234549.6830 

Klaipeda (Lit) 3353590.2730 1302062.9830 5249159.3680 

Kronstadt (Rus) 

1. 1 M/S v It (( id ) 3625340.0140 537853.8050 5202539.0420 

Mantyluoto (Fin) 2831096.7560 1113102.7150 5587165.0170 

Molas (Lit) 3358793.4290 1294907.3540 5247584.3520 

Olands N. U. (Swe) 3295551.6450 1012564.8500 5348113.6740 

Raahe (Fin) 2493889.7970 1131220.2470 5741045.9690 

Ratan (Swe) 2620087.6710 1000008.2270 5709322.5780 

Shepelevo (Rus) 2796394.3740 1556360.1010 5498639.2920 

Spikama (Swe) 2828573.5140 893623.6890 5627447.0490 

Stockholm (Swe) 3101011.4990 1013009.1210 5462375.2820 

Swinoujscie (Pol) 3649458.4600 927709.9170 5130741.6590 

Ustka (Pol) 3545014.3980 1073939.7300 5174949.9340 

Vaasa (Fin) 2691307.2900 1063691.4760 5664806.3350 

Visby (Swe) 3249304.5220 1073624.8500 5364363.1060 

Wamemunde (Ger) 3658217.6960 783004.6440 5148504.2800 

Ger: Germany, Fin: Finland, Swe: Sweden, Lit: Lithuania, Pol: Poland, Rus: Russia_ 


Table 2-3: Cartesian coordinates of GPS stations of BSL third campaign 1997.4, in ITRF 96 reference 
frame. 


Station 

X 

Y 

Z 

Name 

(m) 

(m) 

(m) 

Borkum (Ger) 

3770667.9990 

446076.4900 

5107686.2080 

Degerby (Fin) 

2994064.9360 

1112559.0570 

5502241.3760 

Furuogrund (Swe) 

2527022.8720 

981957.2890 

5753940.9920 

Hamina (Fin) 

2795471.2050 

1435427.7920 

5531682.2000 

Hanko (Fin) 

2959210.9710 

1254679.1200 

5490594.4410 

Helgoland (Ger) 

3706044.9440 

513713.2150 

5148193.4470 

Helsinki (Fin) 

2885137.3910 

1342710.2300 

5509039.1190 

Kemi (Fin) 

2397071.5770 

1093330.3130 

5789108.4470 

Klagshamn (Swe) 

3527585.7670 

807513.8950 

5234549.7020 

Klaipeda (Lit) 

3353590.2430 

1302063.0140 

5249159.4120 

Kronstadt (Rus) 

2776311.8190 

1587590.1310 

5499880.1330 

List/Sylt (Ger) 

3625339.9122 

537853.8700 

5202539.0260 

Mantyluoto (Fin) 

2831096.7190 

1113102.7640 

5587165.0460 

Molas (Lit) 

3358793.3810 

1294907.4050 

5247584.4010 
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Tabelle 2-3: (continued) 


Station 

X 

Y 

Z 

Name 

(m) 

(m) 

(m) 

Olands N. U. (Swe) 

3295551.5710 

1012564.9060 

5348113.6690 

Raahe (Fin) 

2494035.0240 

1131370.9940 

5740955.4100 

Ratan (Swe) 

2620087.6290 

1000008.2700 

5709322.6040 

Shepelevo (Rus) 

2796394.9140 

1556363.7830 

5498638.0600 

Spikama (Swe) 

2828573.4640 

893623.7290 

5627447.0690 

Stockholm (Swe) 

3101008.8620 

1013021.0370 

5462373.3830 

Swinoujscie (Pol) 

3648326.5170 

924984.0310 

5132035.2720 

Ustka (Pol) 

3545014.3300 

1073939.7720 

5174949.9470 

Vaasa (Fin) 

2691307.2540 

1063691.5240 

5664806.3800 

Visby (Swe) 

3249304.4370 

1073624.8910 

5364363.0730 

Wamemtinde (Ger) 

3658231.7070 

783518.3220 

5148404.3509 

Ger: Germany, Fin: Finland, Swe: Sweden, Lit: Lithuania, Pol: Poland, Rus: Russia 


Table 2-4: Jacobi ellipsoidal coordinates of the GPS stations of BSL project first campaign 1990.8, with re¬ 
spect to reference ellipsoid of WGD2000 in mean tide system, a = (6378136.701 ± 0.053)m and 
b = (6356751.661 ± 0.052)m. 


Station 

Name 

o 

A 

> 


o 

<t> 

> 

? 5 

u 

m 

Borkum (Ger) 

6 

44 

48.592 

53 

27 

56.290 

6356797.4401 

Degerby (Fin) 

20 

23 

5.799 

59 

56 

54.812 

6356772.7694 

Furuogrund (Swe) 

21 

14 

6.953 

64 

50 

43.972 

6356784.9708 

Hamina (Fin) 

27 

10 

47.061 

60 

28 

56.235 

6356769.0382 

Hanko (Fin) 

22 

58 

38.020 

59 

44 

21.212 

6356773.7771 

Helgoland (Ger) 

7 

53 

30.345 

54 

5 

0.468 

6356795.7627 

Helsinki (Fin) 

24 

57 

23.340 

60 

4 

14.304 

6356776.0350 

Kemi (Fin) 

24 

30 

56.527 

65 

36 

5.357 

6356772.7375 

Klagshamn (Swe) 

12 

53 

37.154 

55 

25 

56.873 

6356790.0414 

Klaipeda (Lit) 

- 

- 

- 

- 

- 

~ 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

- 

- 

- 

List/Sylt (Ger) 

56 

1 

12.319 

38 

47 

11.504 

8305186.1308 

Mantyluoto (Fin) 

21 

27 

47.774 

61 

30 

49.223 

6356773.4710 

Molas (Lit) 

- 

- 

- 

- 

- 

- 

- 

Olands N. U. (Swe) 

17 

4 

46.851 

57 

16 

48.643 

6356783.6250 

Raahe (Fin) 

24 

24 

52.486 

64 

35 

0.055 

6356772.4568 

Ratan (Swe) 

20 

53 

25.243 

63 

54 

56.307 

6356775.1165 

Shepelevo (Rus) 

- 

- 

- 

- 

- 

- 

- 

Spikama (Swe) 

17 

31 

57.907 

62 

17 

3.822 

6356779.3644 

Stockholm (Swe) 

18 

5 

26.484 

59 

14 

16.253 

6356788.2619 

Swinoujscie (Pol) 

14 

15 

45.951 

53 

48 

58.479 

6356790.0989 

Ustka (Pol) 

16 

51 

13.868 

54 

29 

48.357 

6356786.0870 

Vaasa (Fin) 

21 

33 

55.917 

63 

1 

3.025 

6356771.3215 

Visby (Swe) 

18 

17 

3.922 

57 

33 

7.925 

6356779.3580 

Wamemimde (Ger) 

12 

5 

19.600 

54 

5 

11.799 

6356793.1157 
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Table 2-5: Jacobi ellipsoidal coordinates of the GPS stations of BSL project, second campaign 1993.4, with 
respect to reference ellipsoid of WGD2000 in mean tide system, a = (6378136.701 =b 0.053)m 

b = (6356751.661 ± 0.052)m. 


Station Name 

o 

A 

,, 

O 

<£ 

5 9 

u 

m 

Borkum (Ger) 

6 

44 

48.587 

53 

27 

56.294 

6356797.0766 

Degerby (Fin) 

20 

23 

5.802 

59 

56 

54.816 

6356772.8204 

Furuogrund (Swe) 

21 

14 

6.945 

64 

50 

43.974 

6356785.1749 

Hamina (Fin) 

27 

10 

47.066 

60 

28 

56.239 

6356769.0445 

Hanko (Fin) 

22 

58 

38.023 

59 

44 

21.217 

6356773.8537 

Helgoland (Ger) 

7 

53 

30.344 

54 

5 

0.474 

6356795.8647 

Helsinki (Fin) 

24 

57 

23.341 

60 

4 

14.307 

6356776.1370 

Kemi (Fin) 

24 

30 

56.525 

65 

36 

5.358 

6356772.6609 

Klagshamn (Swe) 

12 

53 

37.155 

55 

25 

56.877 

6356790.2264 

Klaipeda (Lit) 

21 

13 

9.013 

55 

39 

54.157 

6356805.2088 

Kronstadt (Rus) 

- 

- 

- 

- 

- 

- 

- 

List/Sylt (Ger) 

8 

26 

19.755 

54 

55 

37.511 

6356797.0447 

Mantyluoto (Fin) 

21 

27 

47.774 

61 

30 

49.228 

6356773.4901 

Molas (Lit) 

21 

4 

58.889 

55 

38 

24.697 

6356781.6159 

Olands N. U. (Swe) 

17 

4 

46.850 

57 

16 

48.646 

6356783.7334 

Raahe (Fin) 

24 

23 

55.998 

64 

34 

25.438 

6356772.1890 

Ratan (Swe) 

20 

53 

25.235 

63 

54 

56.312 

6356775.1357 

Shepelevo (Rus) 

29 

5 

54.760 

59 

52 

59.790 

6356771.8318 

Spikama (Swe) 

17 

31 

57.902 

62 

17 

3.826 

6356779.5047 

Stockholm (Swe) 

18 

5 

26.485 

59 

14 

16.256 

6356788.4469 

Swinoujscie (Pol) 

14 

15 

45.953 

53 

48 

58.483 

6356790.2519 

Ustka (Pol) 

16 

51 

13.872 

54 

29 

48.361 

6356786.1954 

Vaasa (Fin) 

21 

33 

55.911 

63 

1 

3.028 

6356771.4300 

Visby (Swe) 

18 

17 

3.925 

57 

33 

7.929 

6356779.5685 

Wamemunde (Ger) 

12 

4 

52.651 

54 

5 

16.952 

6356811.9505 


Table 2-6: Jacobi ellipsoidal coordinates of the GPS stations of BSL project, third campaign 1997.4, with re¬ 
spect to reference ellipsoid of WGD2000 in mean tide system, a = (6378136.701 db 0.053)m 

b = (6356751.661 ± 0.052)m . 


Station Name 

o 

A 


O 

<t> 

9 

,, 

u 

m 

Borkum (Ger) 

6 

44 

48.591 

53 

27 

56.296 

6356797.0000 

Degerby (Fin) 

20 

23 

4.091 

59 

56 

52.837 

6356773.9748 

Furuogrund (Swe) 

21 

14 

6.949 

64 

50 

43.975 

6356785.1940 

Hamina (Fin) 

27 

10 

47.069 

60 

28 

56.240 

6356769.0318 

Hanko (Fin) 

22 

58 

35.444 

59 

44 

20.373 

6356777.1894 

Helgoland (Ger) 

7 

53 

30.348 

54 

5 

0.476 

6356795.9221 

Helsinki (Fin) 

24 

57 

24.245 

60 

4 

13.965 

6356776.5070 

Kemi (Fin) 

24 

31 

5.674 

65 

36 

7.401 

6356778.4778 

Klagshamn (Swe) 

12 

53 

37.160 

55 

25 

56.878 

6356790.2200 

Klaipeda (Lit) 

21 

13 

9.016 

55 

39 

54.158 

6356805.2343 

Kronstadt (Rus) 

29 

45 

44.638 

59 

54 

20.012 

6356772.6865 

List/Sylt (Ger) 

8 

26 

19.759 

54 

55 

37.513 

6356796.9809 

Mantyluoto (Fin) 

21 

27 

47.778 

61 

30 

49.229 

6356773.5092 
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Tabelle 2-6: (continued) 


Station Name 

o 

A 


O 

4> 

9 9 

u 

m 

Molas (Lit) 

21 

4 

58.893 

55 

38 

24.698 

6356781.6414 

Olands N. U. (Swe) 

17 

4 

46.854 

57 

16 

48.648 

6356783.6952 

Raahe (Fin) 

24 

24 

1.820 

64 

34 

18.495 

6356773.6451 

Ratan (Swe) 

20 

53 

25.239 

63 

54 

56.313 

102253.32707 

Shepelevo (Rus) 

29 

5 

54.951 

59 

52 

59.707 

6356771.9019 

Spikama (Swe) 

17 

31 

57.906 

62 

17 

3.827 

6356779.5047 

Stockholm (Swe) 

18 

5 

27.253 

59 

14 

16.192 

6356787.4200 

Swinoujscie (Pol) 

14 

13 

36.455 

53 

50 

9.408 

6356794.5253 

Ustka (Pol) 

16 

51 

13.875 

54 

29 

48.363 

6356786.1763 

Vaasa (Fin) 

21 

33 

55.915 

63 

1 

3.029 

6356771.4619 

Visby (Swe) 

18 

17 

3.929 

57 

33 

7.931 

6356779.5047 

Wamemimde (Ger) 

12 

5 

20.183 

54 

5 

11.875 

6356801.8730 


Table 2-7: Orthometric height of the GPS stations of BSL project, first campaign 1990.8, derived from dif¬ 
ferent sources (casel-case5). 


Station 

Name 

orthometric 
height (case 1) 
(Af. Poutanen 
etal. , 1999) 

orthometric 
height (case 2) 
(J. Kakkuri , 
1995) 

orthometric 
height (case 3) 

(J. Kakkuri and 
M. Poutanen , 
1997) 

orthometric 
height 
(case 4) 

(M. Poutanen 
etal , 1999) 

orthometric 
height 
(case 5) 

( EGG97) 

Borkum (Ger) 

4.581 

5.454 

- 

5.884 

5.556 

Degerby (Fin) 

1.695 

2.345 

2.428 

2.541 

2.422 

Furuogrund (Swe) 

10.912 

11.597 

11.541 

11.830 

12.180 

Hamina (Fin) 

1.619 

2.374 

2.357 

2.645 

2.458 

Hanko (Fin) 

1.762 

2.434 

2.473 

2.694 

2.626 

Helgoland (Ger) 

4.539 

- 

- 

5.196 

4.997 

Helsinki (Fin) 

6.033 

6.654 

6.663 

6.933 

6.820 

Kemi (Fin) 

1.246 

2.161 

2.248 

2.546 

2.485 

Klagshamn (Swe) 

2.039 

2.551 

2.169 

2.775 

2.759 

Klaipeda (Lit) 

- 

- 

- 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

- 

List/Sylt (Ger) 

4.160 

4.904 

- 

5.302 

5.135 

Mantyluoto (Fin) 

2.467 

3.208 

3.222 

3.451 

- 

Molas (Lit) 

- 

- 

- 

- 

- 

Olands N. U. 

4.118 

4.493 

4.714 

4.961 

4.729 

(Swe) 

Raahe (Fin) 

2.289 

3.127 

3.168 

3.510 

3.411 

Ratan (Swe) 

1.476 

2.361 

2.268 

2.581 

2.633 

Shepelevo (Rus) 

- 

- 

- 

- 

- 

Spikama (Swe) 

1.872 

2.938 

2.603 

2.809 

2.943 

Stockholm (Swe) 

12.865 

13.532 

13.620 

13.729 

13.699 

Swinoujscie (Pol) 

2.312 

2.877 

2.516 

2.947 

3.004 

Ustka (Pol) 

1.535 

- 

1.720 

2.180 

2.379 

Vaasa (Fin) 

1.128 

1.828 

1.772 

2.131 

2.060 

Visby (Swe) 

1.986 

2.278 

2.547 

2.713 

2.533 

Warnemiinde (Ger) 

- 

- 

2.963 

3.307 

3.246 
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Table 2-8: Orthometric height of the GPS stations of BSL project, second campaign 1993.4, derived from 
different sources (casel-case5). 

Station 

Name 

orthometric 
height 
(case 1) 

(M. Poutanen 
etal , 1999) 

orthometric 
height 
(case 2) 

(7. Kakkuri , 
1995) 

orthometric 
height 
(case 3) 

(7. Kakkuri and 
M. Poutanen , 
1997) 

orthometric 
height 
(case 4) 

(M. Poutanen 
etal., 1999) 

or¬ 
thometric 
height 
(case 5) 
c EGG97) 

Borkum (Ger) 

4.578 

4.433 

- 

4.740 

4.874 

Degerby (Fin) 

1.681 

1.708 

1.791 

1.745 

1.906 

Furuogrund (Swe) 

10.936 

11.074 

11.018 

11.121 

11.308 

Hamina (Fin) 

1.624 

1.710 

1.693 

1.818 

1.982 

Hanko (Fin) 

1.769 

1.824 

1.863 

1.926 

2.086 

Helgoland (Ger) 

4.537 

- 

- 

4.562 

4.687 

Helsinki (Fin) 

6.039 

6.069 

6.078 

6.187 

6.349 

Kemi (Fin) 

1.266 

1.355 

1.442 

1.562 

1.742 

Klagshamn (Swe) 

2.039 

2.060 

1.678 

2.152 

2.286 

Klaipeda (Lit) 

28.209 

- 

- 

28.323 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

- 

List/Sylt (Ger) 

4.159 

3.930 

- 

4.199 

4.330 

Mantyluoto (Fin) 

2.484 

2.533 

2.547 

2.608 

- 

Molas (Lit) 

4.577 

- 

- 

4.664 

4.799 

Olands N. U. (Swe) 

4.122 

3.906 

4.127 

4.230 

4.375 

Raahe (Fin) 

2.086 

2.139 

2.180 

2.337 

2.523 

Ratan (Swe) 

1.500 

1.660 

1.567 

1.698 

1.881 

Shepelevo (Rus) 

- 

4.211 

4.262 

4.317 

- 

Spikama (Swe) 

1.893 

2.365 

2.030 

2.063 

2.237 

Stockholm (Swe) 

12.873 

13.020 

13.108 

13.062 

13.219 

Swinoujscie (Pol) 

2.309 

2.367 

2.006 

2.315 

2.439 

Ustka (Pol) 

1.532 

- 

1.187 

1.521 

1.649 

Vaasa (Fin) 

1.149 

1.227 

1.171 

1.353 

1.551 

Visby (Swe) 

1.986 

1.775 

2.044 

2.065 

2.212 

Warnemiinde (Ger) 

21.291 

- 

21.146 

21.367 

21.492 


Table 2-9: Orthometric height of the GPS stations of BSL project, 3 rd campaign 1997.4, derived from differ¬ 
ent sources (casel-case5). 


Station 

Name 

orthometric 
height 
(case 1) 

( M. Poutanen 
etal , 1999) 

orthometric 
height 
(case 2) 

(7. Kakkuri , 
1995) 

orthometric 
height 
(case 3) 

(7. Kakkuri and 
M. Poutanen , 
1997) 

orthometric 
height 
(case 4) 

(M. Poutanen 
etal , 1999) 

orthometric 
height 
(case 5) 

( Kakkuri , 
2000) 

Borkum (Ger) 

4.574 

4.405 

- 

4.712 

4.846 

Degerby (Fin) 

2.825 

2.877 

2.960 

2.914 

3.075 

Furuogrund (Swe) 

10.972 

11.108 

11.052 

11.155 

11.342 

Hamina (Fin) 

1.631 

1.697 

1.680 

1.805 

1.969 

Hanko (Fin) 

5.118 

5.173 

5.212 

5.275 

5.435 

Helgoland (Ger) 

4.531 

- 

- 

4.610 

4.735 

Helsinki (Fin) 

6.420 

6.455 

6.464 

6.573 

6.735 

Kemi (Fin) 

7.092 

7.185 

7.272 

7.392 

7.572 
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Tabelle 2-9: (continued) 


Station 

Name 

orthometric 
height 
(case 1) 

(M. Poutanen 
etal., 1999) 

orthometric 
height 
(case 2) 

(/. Kakkuri , 
1995) 

orthometric 
height 
(case 3) 

(J. Kakkuri and 
M. Poutanen , 
1997) 

orthometric 
height 
(case 4) 

(M. Poutanen 
etal , 1999) 

orthometric 
height 
(case 5) 

( Kakkuri , 
2000) 

Klagshamn (Swe) 

2.038 

2.099 

1.717 

2.191 

2.325 

Klaipeda (Lit) 

28.209 

- 

- 

28.336 

- 

Kronstadt (Rus) 

- 

- 

- 

4.739 

- 

List/Sylt (Ger) 

4.155 

3.916 

- 

4.185 

4.316 

Mantyluoto (Fin) 

2.509 

2.562 

2.576 

2.637 

- 

Molas (Lit) 

4.577 

- 

- 

4.679 

4.814 

Olands N. U. 

4.127 

3.917 

4.138 

4.241 

4.386 

(Swe) 

Raahe (Fin) 

3.528 

3.607 

3.648 

3.805 

3.991 

Ratan (Swe) 

1.535 

1.691 

1.598 

1.729 

1.912 

Shepelevo (Rus) 

- 

4.282 

4.333 

4.388 

- 

Spikama (Swe) 

1.924 

2.390 

2.055 

2.088 

2.262 

Stockholm (Swe) 

11.905 

12.027 

12.115 

12.069 

12.226 

Swinoujscie (Pol) 

- 

6.681 

6.320 

6.629 

6.753 

Ustka (Pol) 

1.528 

- 

1.181 

1.515 

1.643 

Vaasa (Fin) 

1.180 

1.275 

1.219 

1.401 

1.599 

Visby (Swe) 

1.992 

1.771 

2.040 

2.061 

2.208 

Wamemimde Ger) 

11.319 

- 

11.089 

11.310 

11.435 


3. Baltic Sea Level Project 
First Campaign 1990.8 

Here we present the results of w 0 computation 
based on the GPS observations of the Baltic Sea 
Level project, 1 st campaign, in five cases introduced 
before. Table 3-1-Table 3-5 are presenting the 


gravitational potential U(\,</),u), centrifugal po¬ 
tential V(X ,</>,?/), free-air gravity potential 
A W{\ 4>, u) computed based on the orthometric 

heights of case 1-case 5. A summary of the mean 
values of w Q derived from different cases and the 
standard deviation of those mean values is given in 
Table 3-6. 


Table 3-1: Gravitational potential </>, u ), centrifugal potential V(\<j),u ), free-air gravity potential 
A W( A, 0, u) computed based on the orthometric heights of the case 1 (M. Poutanen et al ., 1999), and gauge 
value of geoid potential W 0 for the epoch of 1990.8. The W 0 at Lisk (Ger) found as outlier, not included in 
the mean. 


Station 

Name 

Gravitational 

Potential 

U (A, <f>, u) 

(m 2 s~ 2 ) 

Centrifugal 

Potential 

V(\</>,u) 

( m 2 s ~ 2 ) 

Free Air 
Reduction 

AW(\,<p,u) 

(m 2 s~ 2 ) 

Gauge 

Value 

Wo 

(m 2 s~ 2 ) 

Borkum (Ger) 

62598475.745 

38330.956 

44.956 

62636851.657 

Degerby (Fin) 

62609712.120 

27124.194 

16.643 

62636852.957 

Furuogrund (Swe) 

62617207.876 

19542.016 

107.183 

62636857.076 

Hamina (Fin) 

62610584.084 

26255.428 

15.897 

62636855.409 

Hanko (Fin) 

62609368.721 

27467.501 

17.301 

62636853.523 
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Tabelle 3-1: (continued) 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

U(\,<t>,u) 

V(\,<f>,u) 

u) 



( m 2 s ~ 2 ) 

(: m 2 s ~ 2 ) 

(■ m 2 s ~ 2 ) 

(mV 2 ) 

Helgoland (Ger) 

62599594.591 

37218.925 

44.546 

62636858.063 

Helsinki (Fin) 

62609873.004 

26924.677 

59.239 

62636856.920 

Kemi (Fin) 

62618386.248 

18455.996 

12.240 

62636854.484 

Klagshamn (Swe) 

62602016.194 

34818.793 

20.014 

62636855.001 

Klaipeda (Lit) 

- 

- 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

47933936.649 

111867.236 

23.866 

48045827.752 

Mantyluoto (Fin) 

62612227.796 

24604.331 

24.226 

62636856.353 

Molas (Lit) 

- 

- 

- 

- 

Olands N. U. (Swe) 

62605215.375 

31601.731 

40.426 

62636857.531 

Raahe (Fin) 

62616906.611 

19924.257 

22.483 

62636853.351 

Ratan (Swe) 

62615931.246 

20910.717 

14.497 

62636856.460 

Shepelevo (Rus) 

- 

- 

- 

- 

Spikama (Swe) 

62613442.333 

23395.357 

18.384 

62636856.074 

Stockholm (Swe) 

62608433.790 

28295.608 

126.314 

62636855.712 

Swinoujscie (Pol) 

62599134.600 

37698.887 

22.690 

62636856.177 

Ustka (Pol) 

62600357.058 

36479.270 

15.066 

62636851.393 

Vaasa (Fin) 

62614579.338 

22265.876 

11.078 

62636856.292 

Visby (Swe) 

62605700.835 

31135.657 

19.497 

62636855.989 

Wamemlinde (Ger) 

62599616.182 

37213.250 

- 

- 




mean ( m 2 / s 2 ) 

62636855.285 




std(m 2 / s 2 ) 

0.4268 


Table 3-2: Gravitational potential centrifugal potential V(X,(/>,u) 9 free-air gravity potential 

AW(X,(/),u) computed based on the orthometric heights of the case 2 (J. Kakkuri , 1995), and gauge value of 
geoid potential W Q for the epoch of 1990.8. The W () at Lisk (Ger) found as outlier, and not included in the mean. 


Station 

Name 

Gravitational 

Potential 

U(\,<j>,u) 

(■ m 2 s ~ 2 ) 

Centrifugal 

Potential 

V (A, <j>, u) 

(m 2 s~ 2 ) 

Free Air 
Reduction 

A VF(A, <f>, u) 

(: m 2 s~ 2 ) 

Gauge Value 

w 0 

0 m 2 s ~ 2 ) 

Borkum (Ger) 

62598475.745 

38330.956 

47.527 

62636854.228 

Degerby (Fin) 

62609712.120 

27124.194 

16.436 

62636852.751 

Furuogrund (Swe) 

62617207.876 

19542.016 

107.999 

62636857.891 

Hamina (Fin) 

62610584.084 

26255.428 

16.595 

62636856.106 

Hanko (Fin) 

62609368.721 

27467.501 

17.625 

62636853.847 

Helgoland (Ger) 

62599594.591 

37218.925 

- 

- 

Helsinki (Fin) 

62609873.004 

26924.677 

59.082 

62636856.763 

Kemi (Fin) 

62618386.248 

18455.996 

12.780 

62636855.024 

Klagshamn (Swe) 

62602016.194 

34818.793 

19.572 

62636854.559 

Klaipeda (Lit) 

- 

- 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

47933936.649 

111867.236 

22.926 

48045826.811 

Mantyluoto (Fin) 

62612227.796 

24604.331 

24.746 

62636856.873 

Molas (Lit) 

- 

- 

- 

- 

Olands N. U. (Swe) 

62605215.375 

31601.731 

37.795 

62636854.900 



National Height Datum, the Gauss-Listing Geoid Level Value 31 


Tabelle 3-2: (continued) 


Station 

Name 

Gravitational 

Potential 

U (A, <j>, u) 

(i mV 2 ) 

Centrifugal 

Potential 

V(X,<p,u) 

(mV 2 ) 

Free Air 
Reduction 

AW(X, <f>, u) 

(mV 2 ) 

Gauge Value 

w Q 

(mV 2 ) 

Raahe (Fin) 

62616906.611 

19924.257 

23.858 

62636854.726 

Ratan (Swe) 

62615931.246 

20910.717 

17.257 

62636859.220 

Shepelevo (Rus) 

- 

- 

- 

- 

Spikama (Swe) 

62613442.333 

23395.357 

22.784 

62636860.473 

Stockholm (Swe) 

62608433.790 

28295.608 

126.608 

62636856.006 

Swinoujscie (Pol) 

62599134.600 

37698.887 

22.406 

62636855.893 

Ustka (Pol) 

62600357.058 

36479.270 

- 

- 

Vaasa (Fin) 

62614579.338 

22265.876 

11.500 

62636856.714 

Visby (Swe) 

62605700.835 

31135.657 

16.326 

62636852.818 

Wamemunde (Ger) 

62599616.182 

37213.250 

- 

- 


mean (m 2 / s 2 ) 62636855.811 

std(ra 2 /s 2 ) 0.5056 


Table 3-3: Gravitational potential U(X,</),u), centrifugal potential V(X free-air gravity potential 
AW(X,<fi,u) computed based on orthometric heights of the case 3 (/. Kakkuri and M. Poutanen , 1997), and 
gauge value of geoid potential W 0 for the epoch of 1990.8. The W 0 at Ustka (Pol) found as outlier, and not 
included in the mean. 


Station 

Name 

Gravitational 

Potential 

U (A, 4>, u) 

(:mV 2 ) 

Centrifugal 

Potential 

V(X,<f>,u) 

(:mV 2 ) 

Free Air 
Reduction 

AW(X, <f>, u) 

(■mV 2 ) 

Gauge Value 

w 0 

(:mV 2 ) 

Borkum (Ger) 

62598475.745 

38330.956 

- 

- 

Degerby (Fin) 

62609712.120 

27124.194 

17.251 

62636853.566 

Furuogrund (Swe) 

62617207.876 

19542.016 

107.449 

62636857.341 

Hamina (Fin) 

62610584.084 

26255.428 

16.428 

62636855.939 

Hanko (Fin) 

62609368.721 

27467.501 

18.008 

62636854.230 

Helgoland (Ger) 

62599594.591 

37218.925 

- 

- 

Helsinki (Fin) 

62609873.004 

26924.677 

59.170 

62636856.852 

Kemi (Fin) 

62618386.248 

18455.996 

13.634 

62636855.879 

Klagshamn (Swe) 

62602016.194 

34818.793 

15.822 

62636850.809 

Klaipeda (Lit) 

- 

- 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

47933936.649 

111867.236 

- 

- 

Mantyluoto (Fin) 

62612227.796 

24604.331 

24.884 

62636857.010 

Molas (Lit) 

- 

- 

- 

- 

Olands N. U. (Swe) 

62605215.375 

31601.731 

39.964 

62636857.070 

Raahe (Fin) 

62616906.611 

19924.257 

24.261 

62636855.129 

Ratan (Swe) 

62615931.246 

20910.717 

16.344 

62636858.307 

Shepelevo (Rus) 

- 

- 

- 

- 

Spikama (Swe) 

62613442.333 

23395.357 

19.494 

62636857.184 

Stockholm (Swe) 

62608433.790 

28295.608 

127.472 

62636856.870 

Swinoujscie (Pol) 

62599134.600 

37698.887 

18.863 

62636852.350 

Ustka (Pol) 

62600357.058 

36479.270 

11.061 

62636847.389 
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Tabelle 3-3: (continued) 


Station 

Name 

Gravitational 

Potential 

U(\4>,u) 

(mV 2 ) 

Centrifugal 
Potential 
V(X ,4>,u) 

(:mV 2 ) 

Free Air 
Reduction 
AW(\,fau) 

( m 2 s ~ 2 ) 

Gauge Value 

Wo 

(mV 2 ) 

Vaasa (Fin) 

62614579.338 

22265.876 

10.950 

62636856.164 

Visby (Swe) 

62605700.835 

31135.657 

18.967 

62636855.459 

Wamemunde (Ger) 

62599616.182 

37213.250 

22.367 

62636851.798 



mean (m 2 / s 2 ) 

62636855.409 




std(m 2 /s 2 ) 

0.5227 


Table 3-4: Gravitational potential centrifugal potential V(A, 0, u ), free-air gravity potential 

AW(\,<f),u) computed based on orthometric heights of the case 4 (Af. Poutanen et al 1999), and gauge 
value of geoid potential W 0 for the epoch of 1990.8. The W 0 at Ustka (Pol) and Lisk (Ger) are outliers and 
excluded from mean. 


Station 

Name 

Gravitational 

Potential 

U{\cj>,u) 

(imV 2 ) 

Centrifugal 

Potential 

V{\cj>,u) 

(mV 2 ) 

Free Air 
Reduction 

(mV 2 ) 

Gauge Value 

w 0 

(m 2 s- 2 ) 

Borkum (Ger) 

62598475.745 

38330.956 

50.540 

62636857.241 

Degerby (Fin) 

62609712.120 

27124.194 

16.800 

62636853.114 

Furubgrund (Swe) 

62617207.876 

19542.016 

108.460 

62636858.352 

Hamina (Fin) 

62610584.084 

26255.428 

17.655 

62636857.167 

Hanko (Fin) 

62609368.721 

27467.501 

18.627 

62636854.849 

Helgoland (Ger) 

62599594.591 

37218.925 

44.782 

62636858.298 

Helsinki (Fin) 

62609873.004 

26924.677 

60.241 

62636857.922 

Kemi (Fin) 

62618386.248 

18455.996 

14.813 

62636857.058 

Klagshamn (Swe) 

62602016.194 

34818.793 

20.475 

62636855.462 

Klaipeda (Lit) 

- 

- 

- 


Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

47933936.649 

111867.236 

24.469 

48045828.354 

Mantyluoto (Fin) 

62612227.796 

24604.331 

25.483 

62636857.609 

Molas (Lit) 

- 

- 

- 

- 

Olands N. U. (Swe) 

62605215.375 

31601.731 

40.975 

62636858.081 

Raahe (Fin) 

62616906.611 

19924.257 

25.803 

62636856.671 

Ratan (Swe) 

62615931.246 

20910.717 

17.631 

62636859.594 

Shepelevo (Rus) 

- 

- 

- 

- 

Spikama (Swe) 

62613442.333 

23395.357 

19.818 

62636857.508 

Stockholm (Swe) 

62608433.790 

28295.608 

127.021 

62636856.419 

Swinoujscie (Pol) 

62599134.600 

37698.887 

21.895 

62636855.382 

Ustka (Pol) 

62600357.058 

36479.270 

14.339 

62636850.667 

Vaasa (Fin) 

62614579.338 

22265.876 

12.738 

62636857.951 

Visby (Swe) 

62605700.835 

31135.657 

19.173 

62636855.665 

Wamemunde (Ger) 

62599616.182 

37213.250 

24.536 
mean(m 2 / s 2 ) 

std(m 2 / s 2 ) 

62636853.967 

62636856.753 

0.3780 
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Table 3-5: Gravitational potential U(A, </>,?/), centrifugal potential V(X,(/>,u) 9 free-air gravity potential 
A TV(A, </>,?/) based on orthometric heights of the case 5 (/. Kakkuri , 2000), and gauge value of geoid poten¬ 
tial W 0 for the epoch of 1990.8. The W 0 at Ustka (Pol) and Lisk (Ger) are outliers and not included in the 
mean value computation. 


Station 

Name 

Gravitational 

Potential 

(mV 2 ) 

Centrifugal 

Potential 

V(\,<f>,u) 

(:mV 2 ) 

Free Air 
Reduction 

(:mV 2 ) 

Gauge 

Value 

(m 2 s~ 2 ) 

Borkum (Ger) 

62598475.745 

38330.956 

51.855 

62636858.556 

Degerby (Fin) 

62609712.120 

27124.194 

18.380 

62636854.695 

Furuogrund (Swe) 

62617207.876 

19542.016 

110.297 

62636860.189 

Hamina (Fin) 

62610584.084 

26255.428 

19.265 

62636858.777 

Hanko (Fin) 

62609368.721 

27467.501 

20.198 

62636856.420 

Helgoland (Ger) 

62599594.591 

37218.925 

46.008 

62636859.525 

Helsinki (Fin) 

62609873.004 

26924.677 

61.831 

62636859.513 

Kemi (Fin) 

62618386.248 

18455.996 

16.581 

62636858.826 

Klagshamn (Swe) 

62602016.194 

34818.793 

21.790 

62636856.777 

Klaipeda (Lit) 

- 

- 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

47933936.649 

111867.236 

25.220 

48045829.106 

Mantyluoto (Fin) 

62612227.796 

24604.331 

- 

- 

Molas (Lit) 

- 

- 

- 

- 

Olands N. U. (Swe) 

62605215.375 

31601.731 

42.399 

62636859.504 

Raahe (Fin) 

62616906.611 

19924.257 

27.630 

62636858.498 

Ratan (Swe) 

62615931.246 

20910.717 

19.428 

62636861.391 

Shepelevo (Rus) 

- 

- 

- 

- 

Spikama (Swe) 

62613442.333 

23395.357 

21.527 

62636859.216 

Stockholm (Swe) 

62608433.790 

28295.608 

128.562 

62636857.960 

Swinoujscie (Pol) 

62599134.600 

37698.887 

23.112 

62636856.599 

Ustka (Pol) 

62600357.058 

36479.270 

15.596 

62636851.923 

Vaasa (Fin) 

62614579.338 

22265.876 

14.682 

62636859.896 

Visby (Swe) 

62605700.835 

31135.657 

20.616 

62636857.108 

Wamemunde (Ger) 

62599616.182 

37213.250 

25.763 
mean(m 2 / s 2 ) 
std(m 2 / s 2 ) 

62636855.194 

62636858.258 

0.4218 


Now let us summarise those results we obtained 
in five different cases for the geoid potential value 
w 0 , based on GPS observations of Baltic Sea level 
project, 1 st campaign, epoch 1990.8, in Table 3-6. 
From a review of Table 3-6 following conclusions 
can be made: 

(i) The most consistent results are 
corresponding to case 4. 

(ii) The results of case 5 show a shift of about 
2.4 (m 2 / s 2 ) with respect to mean value of 
w 0 , calculated from other four cases. 


(Hi) The weighted mean of the geoid potential 
value w 0 based on the GPS observations of 
Baltic Sea Level project, 1 st campaign and 
orthometric heights of all cases is 
(62636856.434 ± 0.558) (m 2 / s 2 ) and 
based on the orthometric heights of casel- 
case4 is 62636855.922 ± 0.366 (m 2 /s 2 ). 
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Table 3-6: Summary of w 0 values computed from the GPS observations of the Baltic Sea level project, first 
campaign 1990.8, for the casel-case5. 


case 

w 0 

(m 2 / s 2 ) 

std(w 0 ) 

(m 2 / s 2 ) 

case 1 

62636855.285 

0.4268 

case 2 

62636855.811 

0.5056 

case 3 

62636855.409 

0.5227 

case 4 

62636856.753 

0.3780 

case 5 

62636858.258 

0.4218 


4. Baltic Sea Level Project, 

Second Campaign 1993.4 

In this section, we will review the results of w Q 
computation based on the GPS observations of Baltic 
Sea level project, 2 nd campaign, in five cases intro¬ 
duced before. Table 4-1-Table 4-5 are presenting 


the gravitational potential centrifugal 

potential V(X,<j),u) 9 free-air gravity potential 
AW(X,c/),u) computed based on the orthometric 

heights of case 1-case 5. A summary of the mean 
values derived from different cases and the standard 
deviation of those mean values is given in Table 4- 
6 . 


Table 4-1: Gravitational potential U(\(/),u), centrifugal potential V(X ,</>, w), free-air gravity potential 
AW(X,(/),u-,H 0 ) , computed based on the orthometric heights of the case 1 (Af. Poutanen et al , 1999), and 
gauge value of geoid potential W 0 for the epoch of 1993.4. 


Station 

Name 

Gravitational 

Potential 

u(\4>,u) 

(: m 2 s - 2 ) 

Centrifugal 

Potential 

V (A, (j>, u) 

(m 2 s~ 2 ) 

Free Air 
Reduction 

SW (A, <j>, u) 

(■mV 2 ) 

Gauge 

Value 

(mV 2 ) 

Borkum (Ger) 

62598479.287 

38330.949 

44.927 

62636855.163 

Degerby (Fin) 

62609711.662 

27124.193 

16.505 

62636852.360 

Furuogrund (Swe) 

62617205.887 

19542.017 

107.419 

62636855.323 

Hamina (Fin) 

62610584.027 

26255.426 

15.946 

62636855.400 

Hanko (Fin) 

62609367.961 

27467.499 

17.370 

62636852.829 

Helgoland (Ger) 

62599593.587 

37218.924 

44.526 

62636857.037 

Helsinki (Fin) 

62609872.036 

26924.677 

59.298 

62636856.010 

Kemi (Fin) 

62618386.991 

18455.995 

12.436 

62636855.423 

Klagshamn (Swe) 

62602014.353 

34818.793 

20.014 

62636853.160 

Klaipeda (Lit) 

62602168.578 

34409.334 

276.885 

62636854.797 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

62601101.969 

35713.293 

40.821 

62636856.083 

Mantyluoto (Fin) 

62612227.594 

24604.328 

24.393 

62636856.316 

Molas (Lit) 

62602356.820 

34452.784 

44.926 

62636854.530 

Olands N. U. (Swe) 

62605214.319 

31601.730 

40.465 

62636856.514 

Raahe (Fin) 

62616895.232 

19938.331 

20.489 

62636854.053 

Ratan (Swe) 

62615931.059 

20910.716 

14.733 

62636856.508 

Shepelevo (Rus) 

62609585.155 

27231.095 

- 

- 

Spikarna (Swe) 

62613440.944 

23395.356 

18.591 

62636854.891 

Stockholm (Swe) 

62608431.963 

28295.609 

126.392 

62636853.964 
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Tabelle 4-1: (continued) 


Station 

Name 

Gravitational 

Potential 

U(X,<j>,u) 

(mV 2 ) 

Centrifugal 

Potential 

v(\ , 0 ,«) 

(mV 2 ) 

Free Air 
Reduction 
8W(\</>,u) 

(■ m 2 s ~ 2 ) 

Gauge 

Value 

(mV 2 ) 

Swinoujscie (Pol) 

62599133.126 

37698.887 

22.661 

62636854.674 

Ustka (Pol) 

62600356.012 

36479.269 

15.036 

62636850.317 

Vaasa (Fin) 

62614578.263 

22265.875 

.284 

62636855.423 

Visby (Swe) 

62605698.773 

31135.657 

19.497 

62636853.927 

Wamemunde (Ger) 

62599433.944 

37210.901 

208.957 

62636853.802 



mean(ra 2 s 2 ) 

62636854.718 




std(ra 2 s 2 ) 

0.3248 


Table 4-2: Gravitational potential U{\,(j),u ), centrifugal potential V(A,</>, u ), free-air gravity potential 
A W(\,(j),u;H 0 ) computed based on orthometric heights of the case 2 (/. Kakkuri , 1995) and gauge value of 
geoid potential W Q for the epoch of 1993.4. 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

U(X, <f>, u) 


SW(X,(j),u) 



(m 2 s ~ 2 ) 

(mV 2 ) 

(■ mV 2 ) 

(mV 2 ) 

Borkum (Ger) 

62598479.287 

38330.949 

43.504 

62636853.740 

Degerby (Fin) 

62609711.662 

27124.193 

16.770 

62636852.625 

Furuogrund (Swe) 

62617205.887 

19542.017 

108.775 

62636856.678 

Hamina (Fin) 

62610584.027 

26255.426 

16.791 

62636856.244 

Hanko (Fin) 

62609367.961 

27467.499 

17.910 

62636853.369 

Helgoland (Ger) 

62599593.587 

37218.924 

- 

- 

Helsinki (Fin) 

62609872.036 

26924.677 

59.592 

62636856.305 

Kemi (Fin) 

62618386.991 

18455.995 

13.310 

62636856.297 

Klagshamn (Swe) 

62602014.353 

34818.793 

20.220 

62636853.366 

Klaipeda (Lit) 

62602168.578 

34409.334 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

62601101.969 

35713.293 

38.574 

62636853.836 

Mantyluoto (Fin) 

62612227.594 

24604.328 

24.874 

62636856.797 

Molas (Lit) 

62602356.820 

34452.784 

- 

- 

Olands N. U. (Swe) 

62605214.319 

31601.730 

38.345 

62636854.393 

Raahe (Fin) 

62616895.232 

19938.331 

21.010 

62636854.573 

Ratan (Swe) 

62615931.059 

20910.716 

16.305 

62636858.079 

Shepelevo (Rus) 

62609585.155 

27231.095 

41.348 

62636857.598 

Spikama (Swe) 

62613440.944 

23395.356 

23.226 

62636859.526 

Stockholm (Swe) 

62608431.963 

28295.609 

127.836 

62636855.407 

Swinoujscie (Pol) 

62599133.126 

37698.887 

23.230 

62636855.243 

Ustka (Pol) 

62600356.012 

36479.269 

- 

- 

Vaasa (Fin) 

62614578.263 

22265.875 

12.050 

62636856.189 

Visby (Swe) 

62605698.773 

31135.657 

17.426 

62636851.855 

Wamemunde (Ger) 

62599433.944 

37210.901 

- 

- 




mean (m 2 s~ 2 ) 

62636855.375 




std (m 2 s~ 2 ) 

0.4557 
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Table 4-3: Gravitational potential U(X,</>,rj) 9 centrifugal potential V(X,<j>,u) 9 free-air gravity potential 
AW(\,(/),u]H 0 ) computed based on orthometric heights of the case 3 (J. Kakkuri and M. Poutanen , 1997), 
and gauge value of geoid potential W 0 for the epoch of 1993.4. The W 0 value at Ustka (Pol) and Klagshamn 
(Swe) are outliers and removed from the mean. 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

U{\Au) 

v{\ 4>, u) 

6W{\, <j>,u) 

w 0 


(mV 2 ) 

(mV 2 ) 

(■ m 2 s - 2 ) 

( m 2 s ~ 2 ) 

Borkum (Ger) 

62598479.287 

38330.949 

- 

- 

Degerby (Fin) 

62609711.662 

27124.193 

17.585 

62636853.440 

Furuogrund (Swe) 

62617205.887 

19542.017 

108.225 

62636856.128 

Hamina (Fin) 

62610584.027 

26255.426 

16.624 

62636856.077 

Hanko (Fin) 

62609367.961 

27467.499 

18.293 

62636853.752 

Helgoland (Ger) 

62599593.587 

37218.924 

- 

- 

Helsinki (Fin) 

62609872.036 

26924.677 

59.681 

62636856.393 

Kemi (Fin) 

62618386.991 

18455.995 

14.165 

62636857.151 

Klagshamn (Swe) 

62602014.353 

34818.793 

16.470 

62636849.616 

Klaipeda (Lit) 

62602168.578 

34409.334 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

62601101.969 

35713.293 

- 

- 

Mantyluoto (Fin) 

62612227.594 

24604.328 

25.012 

62636856.934 

Molas (Lit) 

62602356.820 

34452.784 

- 

- 

Olands N. U. (Swe) 

62605214.319 

31601.730 

40.514 

62636856.563 

Raahe (Fin) 

62616895.232 

19938.331 

21.413 

62636854.976 

Ratan (Swe) 

62615931.059 

20910.716 

15.391 

62636857.166 

Shepelevo (Rus) 

62609585.155 

27231.095 

41.848 

62636858.099 

Spikama (Swe) 

62613440.944 

23395.356 

19.936 

62636856.236 

Stockholm (Swe) 

62608431.963 

28295.609 

128.700 

62636856.271 

Swinoujscie (Pol) 

62599133.126 

37698.887 

19.687 

62636851.700 

Ustka (Pol) 

62600356.012 

36479.269 

11.650 

62636846.931 

Vaasa (Fin) 

62614578.263 

22265.875 

11.500 

62636855.639 

Visby (Swe) 

62605698.773 

31135.657 

20.066 

62636854.496 

Wamemiinde (Ger) 

62599433.944 

37210.901 

207.533 

62636852.379 




mean ( m 2 s ~ 2 ) 

62636855.494 




std(ra 2 s -2 ) 

0.4325 


Table 4-4: Gravitational potential U{X,(j),u) 9 centrifugal potential V(Afree-air gravity potential 
AW(X,(f>,u]H 0 ), based on orthometric heights of the case 4 (Af. Poutanen et al 1999), and gauge value of 
geoid potential W 0 for the epoch of 1993.4. The W 0 at Ustka (Pol) is outlier and removed from the mean. 


Station 

Name 

Gravitational 

Potential 

U(X,4>,u) 

(mV 2 ) 

Centrifugal 

Potential 

V(X ,4>,u) 
(mV 2 ) 

Free Air 
Reduction 
6W(X,(f>,u) 

(m 2 s~ 2 ) 

Gauge 

Value 

W 0 

(: m 2 s- 2 ) 

Borkum (Ger) 

62598479.287 

38330.949 

46.517 

62636856.753 

Degerby (Fin) 

62609711.662 

27124.193 

17.134 

62636852.988 

Furuogrund (Swe) 

62617205.887 

19542.017 

109.236 

62636857.140 

Hamina (Fin) 

62610584.027 

26255.426 

17.851 

62636857.305 

Hanko (Fin) 

62609367.961 

27467.499 

18.911 

62636854.371 
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Tabelle 4-4: (continued) 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

U(\4>,u) 

V(\,<f>,u) 

6W(\,<j>,u) 



(mV 2 ) 

(■ mV 2 ) 

( mV 2 ) 

(mV 2 ) 

Helgoland (Ger) 

62599593.587 

37218.924 

44.772 

62636857.282 

Helsinki (Fin) 

62609872.036 

26924.677 

60.751 

62636857.463 

Kemi (Fin) 

62618386.991 

18455.995 

15.344 

62636858.330 

Klagshamn (Swe) 

62602014.353 

34818.793 

21.123 

62636854.269 

Klaipeda (Lit) 

62602168.578 

34409.334 

278.004 

62636855.916 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

62601101.969 

35713.293 

41.214 

62636856.476 

Mantyluoto (Fin) 

62612227.594 

24604.328 

25.611 

62636857.533 

Molas (Lit) 

62602356.820 

34452.784 

45.779 

62636855.384 

Olands N. U. (Swe) 

62605214.319 

31601.730 

41.525 

62636857.574 

Raahe (Fin) 

62616895.232 

19938.331 

22.955 

62636856.518 

Ratan (Swe) 

62615931.059 

20910.716 

16.678 

62636858.453 

Shepelevo (Rus) 

62609585.155 

27231.095 

42.389 

62636858.639 

Spikama (Swe) 

62613440.944 

23395.356 

20.260 

62636856.560 

Stockholm (Swe) 

62608431.963 

28295.609 

128.248 

62636855.820 

Swinoujscie (Pol) 

62599133.126 

37698.887 

22.720 

62636854.733 

Ustka (Pol) 

62600356.012 

36479.269 

14.928 

62636850.210 

Vaasa (Fin) 

62614578.263 

22265.875 

13.288 

62636857.426 

Visby (Swe) 

62605698.773 

31135.657 

20.273 

62636854.702 

Wamemunde (Ger) 

62599433.944 

37210.901 

209.702 

62636854.548 




mean(ra 2 s~ 2 ) 

62636856.356 




std (m 2 s~ 2 ) 

0.3177 


Table 4-5: Gravitational potential U(\<j),u), centrifugal potential V(\,(j),u ), free-air gravity potential 
AW(\,4>,u;H 0 ), computed based on orthometric heights of the case 5 (/. Kakkuri , 2000), and gauge value 
of geoid potential W 0 for the epoch of 1993.4. The W Q at Ustka (Pol) is outlier and not included in the mean. 


Station 

Name 

Gravitational 

Potential 

U(\,<p,u) 

( m 2 s~ 2 ) 

Centrifugal 

Potential 

V(X ,<f>,u) 

(mV 2 ) 

Free Air 
Reduction 
SW(X,4»,U) 

(m 2 s~ 2 ) 

Gauge 

Value 

w 0 

(■mV 2 ) 

Borkum (Ger) 

62598479.287 

38330.949 

47.832 

62636858.068 

Degerby (Fin) 

62609711.662 

27124.193 

18.714 

62636854.569 

Furubgrund (Swe) 

62617205.887 

19542.017 

111.073 

62636858.977 

Hamina (Fin) 

62610584.027 

26255.426 

19.462 

62636858.915 

Hanko (Fin) 

62609367.961 

27467.499 

20.482 

62636855.942 

Helgoland (Ger) 

62599593.587 

37218.924 

45.999 

62636858.509 

Helsinki (Fin) 

62609872.036 

26924.677 

62.342 

62636859.054 

Kemi (Fin) 

62618386.991 

18455.995 

17.112 

62636860.098 

Klagshamn (Swe) 

62602014.353 

34818.793 

22.438 

62636855.584 

Klaipeda (Lit) 

62602168.578 

34409.334 

- 

- 

Kronstadt (Rus) 

- 

- 

- 

- 

List/Sylt (Ger) 

62601101.969 

35713.293 

42.500 

62636857.762 

Mantyluoto (Fin) 

62612227.594 

24604.328 

- 

- 
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Tabelle 4-5: (continued) 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

U(X,4>,u) 

V (A, <j>, u) 


Wo 


(m 2 s~ 2 ) 

(■mV 2 ) 

(■mV 2 ) 

(mV 2 ) 

Molas (Lit) 

62602356.820 

34452.784 

47.105 

62636856.709 

Olands N. U. (Swe) 

62605214.319 

31601.730 

42.949 

62636858.998 

Raahe (Fin) 

62616895.232 

19938.331 

24.782 

62636858.345 

Ratan (Swe) 

62615931.059 

20910.716 

18.475 

62636860.250 

Shepelevo (Rus) 

62609585.155 

27231.095 

- 

- 

Spikama (Swe) 

62613440.944 

23395.356 

21.969 

62636858.269 

Stockholm (Swe) 

62608431.963 

28295.609 

129.789 

62636857.361 

Swinoujscie (Pol) 

62599133.126 

37698.887 

23.936 

62636855.950 

Ustka (Pol) 

62600356.012 

36479.269 

16.184 

62636851.466 

Vaasa (Fin) 

62614578.263 

22265.875 

15.232 

62636859.371 

Visby (Swe) 

62605698.773 

31135.657 

21.716 

62636856.145 

Wamemunde (Ger) 

62599433.944 

37210.901 

210.929 

62636855.775 




mean (m 2 s~ 2 ) 

62636857.733 




std ( m 2 s ~ 2 ) 

0.3658 


Now let us summarise those results obtained for 
w 0 in five different cases, based on GPS observa¬ 
tion of Baltic Sea level project, 2 nd campaign, epoch 
1993.4, in Table 4-6. From a review of Table 4-6 
following conclusions can be made: 

(i) The results of second campaign are in 
general more accurate than the first 
campaign, which is due to the more accurate 
GPS observation, made under more 
favourable ionspheric condition as is ex¬ 
plained by J. Kakkuri (1995). 

(ii) The most consistent results are corres¬ 
ponding to the case 4. 


(Hi) The results of case 5 show a shift of about 
2.5 (m 2 / s 2 ) with respect to the mean value 
of w 0 calculated from casel-case4. 

(iv) The mean of the geoid potential value w 0 
based on GPS observations of Baltic Sea 
Level project, 2 nd campaign, and the 
orthometric heights of all cases is 
(62636855.962 ± 0.536) (m 2 /s 2 ) and 
based on the orthometric heights of the 
casel-case4 is 62636855.515 =b 
0.385 (m 2 /s 2 ). 


Table 4-6: Summary of w 0 values computed from the GPS observations of the Baltic Sea level project, sec¬ 
ond campaign 1993.4, for the casel-case5. 


case 

w 0 

(m 2 / s 2 ) 

std(w 0 ) 

(m 2 / s 2 ) 

case 1 

62636854.718 

0.3248 

case 2 

62636855.375 

0.4557 

case 3 

62636855.494 

0.4325 

case 4 

62636856.356 

0.3177 

case 5 

62636857.733 

0.3658 



National Height Datum, the Gauss-Listing Geoid Level Value 39 


5. Baltic Sea Level Project, 

Third Campaign 1997.4 

Here we present the results of w 0 computation 
based on the GPS observations of the Baltic Sea 
Level project, 3 rd campaign, in five cases introduced 
before. Table 5-1-Table 5-5 are presenting gravi¬ 


tational potential U{\(j>,u ), centrifugal potential 
V(X,(j>,u) 9 free-air gravity potential A VL(A, 0, u) 

computed based on the orthometric heights of case 
1-case 5. A summary of the mean values derived 
from different cases and the standard deviation of 
those mean values is given in Table 5-6. 


Table 5-1: Gravitational potential U(X,(j>,u), centrifugal potential V(X free-air gravity potential 
A W(X, 0, u;H 0 ) , computed based on orthometric heights of the case 1 (M. Poutanen et al , 1999), and gauge 
value of geoid potential W 0 for the epoch of 1997.4. The W 0 at Ustka (Pol) is outlier and removed from the 
mean. 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

U{\<l>,u) 

V(X 

6W(\, <j>, u) 



(‘ m 2 s ~ 2 ) 

(mV 2 ) 

(■mV 2 ) 

(mV 2 ) 

Borkum (Ger) 

62598480.049 

38330.947 

44.888 

62636855.884 

Degerby (Fin) 

62609699.381 

27125.102 

27.738 

62636852.221 

Furuogrund (Swe) 

62617205.722 

19542.016 

107.773 

62636855.511 

Hamina (Fin) 

62610584.135 

26255.426 

16.015 

62636855.576 

Hanko (Fin) 

62609334.840 

27467.913 

50.254 

62636853.007 

Helgoland (Ger) 

62599593.036 

37218.924 

44.468 

62636856.427 

Helsinki (Fin) 

62609868.216 

26924.835 

63.039 

62636856.090 

Kemi (Fin) 

62618330.677 

18455.223 

69.665 

62636855.565 

Klagshamn (Swe) 

62602014.413 

34818.792 

20.004 

62636853.208 

Klaipeda (Lit) 

62602168.314 

34409.333 

276.885 

62636854.532 

Kronstadt (Rus) 

62609616.854 

27194.595 

- 

- 

List/Sylt (Ger) 

62601102.613 

35713.291 

40.782 

62636856.686 

Mantyluoto (Fin) 

62612227.421 

24604.328 

24.638 

62636856.388 

Molas (Lit) 

62602356.570 

34452.784 

44.926 

62636854.279 

Olands N. U. (Swe) 

62605214.649 

31601.729 

40.514 

62636856.892 

Raahe (Fin) 

62616878.076 

19941.164 

34.653 

62636853.893 

Ratan (Swe) 

62615930.933 

20910.715 

15.077 

62636856.725 

Shepelevo (Rus) 

62609584.472 

27231.134 

- 

- 

Spikama (Swe) 

62613440.933 

23395.355 

18.895 

62636855.183 

Stockholm (Swe) 

62608441.995 

28295.629 

116.888 

62636854.512 

Swinoujscie (Pol) 

62599126.796 

37663.497 

- 

- 

Ustka (Pol) 

62600356.210 

36479.268 

14.997 

62636850.475 

Vaasa (Fin) 

62614577.940 

22265.875 

11.589 

62636855.403 

Visby (Swe) 

62605699.404 

31135.656 

19.556 

62636854.616 

Wamemunde (Ger) 

62599530.263 

37213.313 

111.089 

62636854.665 




mean(ra 2 s~ 2 ) 

62636855.108 




std ( m 2 s ~ 2 ) 

0.3239 
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Table 5-2: Gravitational potential £/(A, </>,?/), centrifugal potential V(A,0, u ), free-air gravity potential 
AW(\,(j),u;H Q ), computed based on orthometric heights of the case 2 (J. Kakkuri , 1995), and gauge value 
of geoid potential W 0 for the epoch of 1997.4. 


Station 

Name 

Gravitational 

Potential 

U(X,<p,u) 

(mV 2 ) 

Centrifugal 

Potential 

u) 

(mV 2 ) 

Free Air 
Reduction 
6W(\<l>,u) 

( m 2 s ~ 2 ) 

Gauge 

Value 

W 0 

(: m 2 s ~ 2 ) 

Borkum (Ger) 

62598480.049 

38330.947 

43.229 

62636854.226 

Degerby (Fin) 

62609699.381 

27125.102 

28.248 

62636852.731 

Furuogrund (Swe) 

62617205.722 

19542.016 

109.109 

62636856.847 

Hamina (Fin) 

62610584.135 

26255.426 

16.663 

62636856.224 

Hanko (Fin) 

62609334.840 

27467.913 

50.794 

62636853.547 

Helgoland (Ger) 

62599593.036 

37218.924 

- 

- 

Helsinki (Fin) 

62609868.216 

26924.835 

63.383 

62636856.433 

Kemi (Fin) 

62618330.677 

18455.223 

70.579 

62636856.478 

Klagshamn (Swe) 

62602014.413 

34818.792 

20.602 

62636853.807 

Klaipeda (Lit) 

62602168.314 

34409.333 

- 

- 

Kronstadt (Rus) 

62609616.854 

27194.595 

- 

- 

List/Sylt (Ger) 

62601102.613 

35713.291 

38.436 

62636854.340 

Mantyluoto (Fin) 

62612227.421 

24604.328 

25.159 

62636856.908 

Molas (Lit) 

62602356.570 

34452.784 

- 

- 

Olands N. U. (Swe) 

62605214.649 

31601.729 

38.453 

62636854.830 

Raahe (Fin) 

62616878.076 

19941.164 

35.429 

62636854.669 

Ratan (Swe) 

62615930.933 

20910.715 

16.609 

62636858.257 

Shepelevo (Rus) 

62609584.472 

27231.134 

42.045 

62636857.650 

Spikama (Swe) 

62613440.933 

23395.355 

23.471 

62636859.760 

Stockholm (Swe) 

62608441.995 

28295.629 

118.086 

62636855.710 

Swinoujscie (Pol) 

62599126.796 

37663.497 

65.568 

62636855.860 

Ustka (Pol) 

62600356.210 

36479.268 

- 

- 

Vaasa (Fin) 

62614577.940 

22265.875 

12.522 

62636856.336 

Visby (Swe) 

62605699.404 

31135.656 

17.386 

62636852.446 

Wamemunde (Ger) 

62599530.263 

37213.313 

- 

- 



mean (m 2 s 2 ) 

62636855.635 




std(ra 2 s 2 ) 

0.4350 


Table 5-3: Gravitational potential U(\(j),u), centrifugal potential V(X free-air gravity potential 
AW(X 1 (j),u;H 0 ), computed based on orthometric heights of the case 3 (/. Kakkuri and M. Poutanen , 1997), 
and gauge value of geoid potential W 0 for the epoch of 1997.4. The W 0 at Ustka (Pol) is outlier and is not 
used for the calculation of mean value. 


Station 

Name 

Gravitational 

Potential 

u{\4>,u) 

(m 2 s~ 2 ) 

Centrifugal 

Potential 

V(\,<j>,u) 

(mV 2 ) 

Free Air 
Reduction 

(mV 2 ) 

Gauge 

Value 

w 0 

(mV 2 ) 

Borkum (Ger) 

62598480.049 

38330.947 

- 

- 

Degerby (Fin) 

62609699.381 

27125.102 

29.063 

62636853.546 

Furuogrund (Swe) 

62617205.722 

19542.016 

108.558 

62636856.296 

Hamina (Fin) 

62610584.135 

26255.426 

16.496 

62636856.057 

Hanko (Fin) 

62609334.840 

27467.913 

51.177 

62636853.930 
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Tabelle 5-3: (continued) 


Station 

Name 

Gravitational 

Potential 

U (A, cf>, u) 

0 m 2 s - 2 ) 

Centrifugal 

Potential 

V{\ <t>, u) 

0 rnh- 2 ) 

Free Air 
Reduction 

SW{\(j>,u) 

(mV 2 ) 

Gauge 

Value 

Wo 

(mV 2 ) 

Helgoland (Ger) 

62599593.036 

37218.924 

- 

- 

Helsinki (Fin) 

62609868.216 

26924.835 

63.471 

62636856.522 

Kemi (Fin) 

62618330.677 

18455.223 

71.433 

62636857.333 

Klagshamn (Swe) 

62602014.413 

34818.792 

16.853 

62636850.058 

Klaipeda (Lit) 

62602168.314 

34409.333 

- 

- 

Kronstadt (Rus) 

62609616.854 

27194.595 

- 

- 

List/Sylt (Ger) 

62601102.613 

35713.291 

- 

- 

Mantyluoto (Fin) 

62612227.421 

24604.328 

25.296 

62636857.046 

Molas (Lit) 

62602356.570 

34452.784 

- 

- 

Olands N. U. (Swe) 

62605214.649 

31601.729 

40.622 

62636857.000 

Raahe (Fin) 

62616878.076 

19941.164 

35.832 

62636855.072 

Ratan (Swe) 

62615930.933 

20910.715 

15.696 

62636857.344 

Shepelevo (Rus) 

62609584.472 

27231.134 

42.546 

62636858.151 

Spikama (Swe) 

62613440.933 

23395.355 

20.182 

62636856.470 

Stockholm (Swe) 

62608441.995 

28295.629 

118.950 

62636856.574 

Swinoujscie (Pol) 

62599126.796 

37663.497 

62.025 

62636852.317 

Ustka (Pol) 

62600356.210 

36479.268 

11.591 

62636847.069 

Vaasa (Fin) 

62614577.940 

22265.875 

11.972 

62636855.786 

Visby (Swe) 

62605699.404 

31135.656 

20.027 

62636855.087 

Wamemiinde (Ger) 

62599530.263 

37213.313 

108.831 

62636852.407 



mean (m 2 s 2 ) 

62636855.389 




std(ra 2 s 2 ) 

0.5059 


Table 5-4: Gravitational potential U(\(j),u), centrifugal potential V(X,0,u), free-air gravity potential 
» based on orthometric heights of the case 4 (M. Poutanen et al. 1999), and gauge value of 
geoid potential W 0 for the epoch of 1997.4. The W Q at Ustka (Pol) is outlier and removed from the mean. 


Station 

Name 

Gravitational 

Potential 

0mV 2 ) 

Centrifugal 

Potential 

v(\ ,<f>,u) 

(m 2 s~ 2 ) 

Free Air 
Reduction 
SW(X,<j>,u) 

(:mV 2 ) 

Gauge 

Value 

W 0 

(mV 2 ) 

Borkum (Ger) 

62598480.049 

38330.947 

46.242 

62636857.238 

Degerby (Fin) 

62609699.381 

27125.102 

28.611 

62636853.095 

Furuogrund (Swe) 

62617205.722 

19542.016 

109.570 

62636857.308 

Hamina (Fin) 

62610584.135 

26255.426 

17.724 

62636857.284 

Hanko (Fin) 

62609334.840 

27467.913 

51.795 

62636854.548 

Helgoland (Ger) 

62599593.036 

37218.924 

45.243 

62636857.203 

Helsinki (Fin) 

62609868.216 

26924.835 

64.541 

62636857.592 

Kemi (Fin) 

62618330.677 

18455.223 

72.612 

62636858.512 

Klagshamn (Swe) 

62602014.413 

34818.792 

21.505 

62636854.710 

Klaipeda (Lit) 

62602168.314 

34409.333 

278.131 

62636855.779 

Kronstadt (Rus) 

62609616.854 

27194.595 

46.533 

62636857.981 

List/Sylt (Ger) 

62601102.613 

35713.291 

41.077 

62636856.980 

Mantyluoto (Fin) 

62612227.421 

24604.328 

25.895 

62636857.645 
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Tabelle 5-4: (continued) 



Gravitational 

Centrifugal 

Free Air 

Gauge 

Station 

Potential 

Potential 

Reduction 

Value 

Name 

u) 

v{\,4>, u) 

6W(X,4>,u) 

w 0 


o m 2 s - 2 ) 

0 m 2 s - 2 ) 

(mV 2 ) 

(■ mV 2 ) 

Molas (Lit) 

62602356.570 

34452.784 

45.927 

62636855.280 

Olands N. U. (Swe) 

62605214.649 

31601.729 

41.633 

62636858.011 

Raahe (Fin) 

62616878.076 

19941.164 

37.374 

62636856.614 

Ratan (Swe) 

62615930.933 

20910.715 

16.982 

62636858.631 

Shepelevo (Rus) 

62609584.472 

27231.134 

43.086 

62636858.691 

Spikama (Swe) 

62613440.933 

23395.355 

20.506 

62636856.794 

Stockholm (Swe) 

62608441.995 

28295.629 

118.498 

62636856.122 

Swinoujscie (Pol) 

62599126.796 

37663.497 

65.057 

62636855.350 

Ustka (Pol) 

62600356.210 

36479.268 

14.869 

62636850.347 

Vaasa (Fin) 

62614577.940 

22265.875 

13.759 

62636857.574 

Visby (Swe) 

62605699.404 

31135.656 

20.233 

62636855.293 

Wamemunde (Ger) 

62599530.263 

37213.313 

111.000 

62636854.576 




mean(ra 2 s~ 2 ) 

62636856.617 




std ( m 2 s ~ 2 ) 

0.2980 


Table 5-5: Gravitational potential U(X,(/),u), centrifugal potential V(X ,0,w), free-air gravity potential 
AW(X,(f),u]H 0 ), based on orthometric heights of the case 5 (/. Kakkuri , 2000), and gauge value of geoid 
potential W 0 for the epoch of 1997.4. The W 0 at Ustka (Pol) is outlier and not included in mean. 


Station 

Name 

Gravitational 

Potential 

(mV 2 ) 

Centrifugal 

Potential 

V(\,<j>,u) 

( m 2 s ~ 2 ) 

Free Air 
Reduction 

(m 2 s~ 2 ) 

Gauge 

Value 

(m 2 s~ 2 ) 

Borkum (Ger) 

62598480.049 

38330.947 

47.557 

62636858.553 

Degerby (Fin) 

62609699.381 

27125.102 

30.192 

62636854.675 

Furudgrund (Swe) 

62617205.722 

19542.016 

111.407 

62636859.145 

Hamina (Fin) 

62610584.135 

26255.426 

19.334 

62636858.895 

Hanko (Fin) 

62609334.840 

27467.913 

53.366 

62636856.119 

Helgoland (Ger) 

62599593.036 

37218.924 

46.470 

62636858.429 

Helsinki (Fin) 

62609868.216 

26924.835 

66.132 

62636859.183 

Kemi (Fin) 

62618330.677 

18455.223 

74.380 

62636860.280 

Klagshamn (Swe) 

62602014.413 

34818.792 

22.821 

62636856.025 

Klaipeda (Lit) 

62602168.314 

34409.333 

- 

- 

Kronstadt (Rus) 

62609616.854 

27194.595 

- 

- 

List/Sylt (Ger) 

62601102.613 

35713.291 

42.362 

62636858.266 

Mantyluoto (Fin) 

62612227.421 

24604.328 

- 

- 

Molas (Lit) 

62602356.570 

34452.784 

47.252 

62636856.605 

Olands N. U. (Swe) 

62605214.649 

31601.729 

43.057 

62636859.434 

Raahe (Fin) 

62616878.076 

19941.164 

39.201 

62636858.441 

Ratan (Swe) 

62615930.933 

20910.715 

18.780 

62636860.428 

Shepelevo (Rus) 

62609584.472 

27231.134 

- 

- 

Spikama (Swe) 

62613440.933 

23395.355 

22.214 

62636858.503 

Stockholm (Swe) 

62608441.995 

28295.629 

120.040 

62636857.664 
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Tabelle 5-5: (continued) 


Station 

Name 

Gravitational 

Potential 

U(X,4>,u) 

0 m 2 s - 2 ) 

Centrifugal 

Potential 

V(\,<f>,u) 

0 mV 2 ) 

Free Air 
Reduction 

AW(X,4>,u-,h 

(mV 2 ) 

Gauge 

Value 

(mV 2 ) 

Swinoujscie (Pol) 

62599126.796 

37663.497 

66.274 

62636856.567 

Ustka (Pol) 

62600356.210 

36479.268 

16.126 

62636851.604 

Vaasa (Fin) 

62614577.940 

22265.875 

15.704 

62636859.518 

Visby (Swe) 

62605699.404 

31135.656 

21.676 

62636856.736 

Wamemiinde (Ger) 

62599530.263 

37213.313 

112.227 

62636855.803 



mean (m 2 s 2 ) 

62636857.963 




std(ra 2 s 2 ) 

0.3590 


Now let us summarise those results we obtained 
in five different cases for w Q , based on GPS obser¬ 
vation of the Baltic Sea Level project, 3 rd campaign, 
epoch 1997.4, in Table 5-6. From a review of Ta¬ 
ble 5-6 following conclusions can be made: 

(i) The results of w 0 computations based on the 
data of the third campaign are in general 
more accurate than the results we obtained 
for the first and the second campaigns, 
which indicates the improvement in the 
accuracy of the GPS observations. 

(ii) The most consistent results are 
corresponding to the heights of case 4. 

(iii) The results of case 5 shows a shift of about 
2.3 (m 2 / s 2 ) with respect to mean value of 
w 0 calculated from other four cases. 


The weighted mean of the geoid potential 
value w 0 based on the GPS observations of 
Baltic Sea Level project, 3 rd campaign, and 
orthometric heights of all cases is 
(62636856.248 ± 0.525) (m 2 / s 2 ), and 
based on the orthometric heights of the 
casel-case4 is 62636855.804 ± 0.378 
(m 2 /* 2 ). 

Finally, from a review of Table 3-6 , Table 
4-6 , and Table 5-6 we can conclude the 
value w 0 = (62636855.75 ±0.21)(m 2 / s 2 ) 

as our best estimate of the geoid potential 
value. 


Table 5-6: Summary of w 0 values computed from the GPS observations of the Baltic Sea level project, third 
campaign 1997.4, for casel-case5. 


Case 

w 0 

( m 2 / s 2 ) 

std(w Q ) 

(m 2 / s 2 ) 

case 1 

62636855.108 

0.3239 

case 2 

62636855.635 

0.4350 

case 3 

62636855.389 

0.5059 

case 4 

62636856.617 

0.2980 

case 5 

62636857.963 

0.3590 
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6. Time Derivative of W 0 

In the Previous section we presented the value 
w 0 = (62636855.75 ±0.21)(m 2 / s 1 ), from the 

analysis of w 0 values computed in five cases and 
three campaigns of GPS observations of Baltic Sea 
Level project. However, for w Q computation we 
need a different computation strategy. This strategy 
is motivated by the fact that in relatively short time 
span from 1990.8 till 1997.4, the only variation in 
w 0 which can be sensible, is due to eustatic rise. 
Therefore, to remove the effects of any sources 
other than variations in sea surface, we computed 
again the w Q values for the epochs of 1990.8 and 
1993.4 using the orthometric heights of casel- 
case5 but the GPS coordinates of Baltic Sea Level 
project, 3 rd campaign, epoch 1997.4 only, which are 
the most accurate GPS observations being made so 
far at the Baltic Sea. 


Based on the new computation, we have come 

up with the results shown in Table 6-1-Table 6-5. 

From a review of Table 6-1-Table 6-5 following 

conclusions can be made: 

(i) The most consistent results correspond to 
case 1, where the heights are computed in 
their respective national heights without 
application of regional geoid solutions. 

(ii) Within the results of case 1, most consistent 
results belong to the tide gauge stations of 
Germany. 

(iii) From the w Q computations in four stations 
Borkum (Ger), Helgoland (Ger), List/Sylt 
(Ger), Ustka (Pol) we have computed 
following estimatation for the time 
derivative of the geoid potential value “ w 0 ”. 

w 0 - (-0.0099 ± 0.00079)(m 2 /s 2 )/year (6.1) 


Table 6-1: Results of w 0 computation based on the orthometric heights of case 1. The most consistent re¬ 
sults correspond to tide gauge stations Borkum (Ger), Helgoland (Ger), List/Sylt (Ger), Ustka (Pol). 


Station 

Name 

w 0 

(m 2 s~ 2 )/ y 

1993.4-1990.8 

w 0 

(mV 2 )/ y 

1997.4-1993.4 

w 0 

(m 2 s~ 2 )/y 

1997.4-1990.8 

w 0 

(mV 2 )/ y 

mean 

Borkum(Ger) 

-0.01153 

-0.00975 

-0.01045 

-0.01058 

Degerby(Fin) 

-0.05307 

2.80830 

1.68110 

1.47870 

Furuogrund(Swe) 

0.09076 

0.08850 

0.08939 

0.08955 

Hamina(Fin) 

0.01884 

0.01725 

0.01787 

0.01799 

Hanko(Fin) 

0.02653 

8.22100 

4.99290 

4.41350 

Helgoland(Ger) 

-0.00769 

-0.01475 

-0.01197 

-0.01147 

Helsinki(Fin) 

0.02269 

0.93525 

0.57576 

0.51123 

Kemi(Fin) 

0.07576 

14.30700 

8.70090 

7.69460 

Klagshamn(Swe) 

0 

-0.00250 

-0.00151 

-0.00133 

Klaipeda(Lit) 

- 

0 

- 

- 

Kronstadt(Rus) 

- 

- 

- 

- 

List/Sylt(Ger) 

-0.00384 

-0.00975 

-0.00742 

-0.00700 

Mantyluoto(Fin) 

0.06423 

0.06150 

0.06257 

0.06276 

Molas(Lit) 

- 

0 

- 

- 

OlandsN.U.(Swe) 

0.01538 

0.01225 

0.01348 

0.01370 

Raahe(Fin) 

-0.76692 

3.54100 

1.84390 

1.53930 

Ratan(Swe) 

0.09038 

0.08600 

0.08772 

0.08803 

Shepelevo(Rus) 

- 

- 

- 

- 

Spikama(Swe) 

0.07923 

0.07600 

0.07727 

0.07750 

Stockholm(Swe) 

0.03000 

-2.37600 

-1.42820 

-1.25810 

Swinoujscie(Pol) 

-0.01115 

- 

- 

- 

Ustka(Pol) 

-0.01153 

-0.00975 

-0.01045 

-0.01058 

Vaasa(Fin) 

0.07923 

0.07600 

0.07727 

0.07750 

Visby(Swe) 

0 

0.01475 

0.00893 

0.00789 

W amemimde(Ger) 

- 

-24.46700 

- 

- 
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Table 6-2: Results of w 0 computation based on the orthometric heights of case 2. The results from one sta¬ 
tion to anther are very much different and therefore not acceptable. 


Station 

Name 

w 0 

('m 2 s~ 2 )/y 

1993.4-1990.8 

w 0 

(mV 2 )/*/ 

1997.4-1993.4 

w 0 

(mV 2 )/*/ 

1997.4-1990.8 

w 0 

(mV 2 )/*/ 

mean 

Borkum(Ger) 

-1.5477 

-0.0685 

-0.65121 

-0.7558 

Degerby(Fin) 

0.12846 

2.8695 

1.7897 

1.5959 

Furuogrund(Swe) 

0.29846 

0.0835 

0.16818 

0.18338 

Hamina(Fin) 

0.075385 

-0.03175 

0.010455 

0.01803 

Hanko(Fin) 

0.10962 

8.221 

5.0256 

4.4521 

Helgoland(Ger) 

- 

- 

- 


Helsinki(Fin) 

0.19615 

0.9475 

0.65152 

0.59839 

Kemi(Fin) 

0.20385 

14.317 

8.7573 

7.7594 

Klagshamn(Swe) 

0.24923 

0.09575 

0.15621 

0.16706 

Klaipeda(Lit) 

- 


- 

- 

Kronstadt(Rus) 

- 

- 

- 

- 

List/Sylt(Ger) 

-0.24923 

-0.03425 

-0.11894 

-0.13414 

Mantyluoto(Fin) 

0.048846 

0.07125 

0.062424 

0.06084 

Molas(Lit) 

- 

- 

- 

- 

OlandsN.U.(Swe) 

0.21154 

0.027 

0.099697 

0.11275 

Raahe(Fin) 

-1.0958 

3.6047 

1.753 

1.4207 

Ratan(Swe) 

-0.36654 

0.076 

-0.098333 

-0.12962 

Shepelevo(Rus) 

- 

0.17425 

- 

- 

Spikama(Swe) 

0.17 

0.0615 

0.10424 

0.11191 

Stockholm(Swe) 

0.47231 

-2.4375 

-1.2912 

-1.0855 

Swinoujscie(Pol) 

0.31692 

10.584 

6.5397 

5.8137 

Ustka(Pol) 

- 

- 

- 

- 

Vaasa(Fin) 

0.21154 

0.11775 

0.1547 

0.16133 

Visby(Swe) 

0.42269 

-0.00975 

0.16061 

0.19118 

W amemunde(Ger) 

- 

- 

- 

- 


Table 6-3: Results of w Q computation based on the orthometric heights of case 3. The results from one sta¬ 
tion to anther are very much different and therefore not acceptable. 


Station 

Name 

w 0 

(mV 2 )/*/ 

1993.4-1990.8 

w 0 

(mV 2 )/ y 

1997.4-1993.4 

w 0 

(' m 2 s~ 2 )/y 

1997.4-1990.8 

w 0 

(m 2 s~ 2 ) / y 

mean 

Borkum(Ger) 

- 

- 

- 

- 

Degerby(Fin) 

0.12846 

2.8695 

1.7897 

1.5959 

Furuogrund(Swe) 

0.29846 

0.08325 

0.16803 

0.18325 

Hamina(Fin) 

0.075769 

-0.032 

0.010455 

0.018075 

Hanko(Fin) 

0.10962 

8.221 

5.0256 

4.4521 

Helgoland(Ger) 

- 

- 

- 

- 

Helsinki(Fin) 

0.19615 

0.94775 

0.65167 

0.59852 

Kemi(Fin) 

0.20423 

14.317 

8.7574 

7.7596 

Klagshamn(Swe) 

0.24923 

0.09575 

0.15621 

0.16706 

Klaipeda(Lit) 

- 

- 

- 

- 

Kronstadt(Rus) 

- 

- 

- 

- 

List/Sylt(Ger) 

- 

- 

- 

- 
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Tabelle 6-3: (continued) 


Station 

Name 

w 0 

(imV 2 )/?/ 

1993.4-1990.8 

(mV 2 )/ y 

1997.4-1993.4 

w 0 

(■ m 2 s~ 2 )/y 

1997.4-1990.8 

w 0 

(m 2 s~ 2 )/y 

mean 

Mantyluoto(Fin) 

0.049231 

0.07125 

0.062576 

0.061019 

Molas(Lit) 

- 

- 

- 

- 

OlandsN.U.(Swe) 

0.21154 

0.027 

0.099697 

0.11275 

Raahe(Fin) 

-1.0954 

3.6047 

1.7532 

1.4208 

Ratan(Swe) 

-0.36654 

0.07625 

-0.098182 

-0.12949 

Shepelevo(Rus) 

- 

0.17425 

- 

- 

Spikama(Swe) 

0.17 

0.0615 

0.10424 

0.11191 

Stockholm(Swe) 

0.47231 

-2.4375 

-1.2912 

-1.0855 

Swinoujscie(Pol) 

0.31692 

10.584 

6.5397 

5.8137 

Ustka(Pol) 

0.22654 

-0.01475 

0.080303 

0.097364 

Vaasa(Fin) 

0.21154 

0.11775 

0.1547 

0.16133 

Visby(Swe) 

0.42269 

-0.00975 

0.16061 

0.19118 

W amemiinde(Ger) 

71.218 

-24.676 

13.101 

19.881 


Table 6-4: Results of w 0 computation based on the orthometric heights of case 4. The results from one sta¬ 
tion to anther are very much different and therefore not acceptable. 


Station 

Name 

w 0 

(m 2 s~ 2 )/ y 

1993.4-1990.8 

w 0 

(m 2 s~ 2 )/ y 

1997.4-1993.4 

w 0 

(m 2 s~ 2 )/ y 

1997.4-1990.8 

w 0 

( m 2 s~ 2 )/y 

mean 

Borkum(Ger) 

-1.5477 

-0.06875 

-0.65136 

-0.75594 

Degerby(Fin) 

0.12846 

2.8695 

1.7897 

1.5959 

Furuogrund(Swe) 

0.29846 

0.0835 

0.16818 

0.18338 

Hamina(Fin) 

0.075385 

-0.032 

0.010303 

0.017896 

Hanko(Fin) 

0.10962 

8.2208 

5.0255 

4.4519 

Helgoland(Ger) 

-0.0038462 

0.118 

0.07 

0.061385 

Helsinki(Fin) 

0.19654 

0.9475 

0.65167 

0.59857 

Kemi(Fin) 

0.20385 

14.317 

8.7574 

7.7595 

Klagshamn(Swe) 

0.24923 

0.09575 

0.15621 

0.16706 

Klaipeda(Lit) 

- 

0.032 

- 

- 

Kronstadt(Rus) 

- 

- 

- 

- 

ListZSylt(Ger) 

-0.24885 

-0.0345 

-0.11894 

-0.1341 

Mantyluoto(Fin) 

0.049231 

0.07125 

0.062576 

0.061019 

Molas(Lit) 

- 

0.03675 

- 

- 

OlandsN.U.(Swe) 

0.21154 

0.027 

0.099697 

0.11275 

Raahe(Fin) 

-1.0954 

3.6047 

1.7532 

1.4208 

Ratan(Swe) 

-0.36654 

0.07625 

-0.098182 

-0.12949 

Shepelevo(Rus) 

- 

0.17425 

- 

- 

Spikama(Swe) 

0.16962 

0.0615 

0.10409 

0.11174 

Stockholm(Swe) 

0.47192 

-2.4375 

-1.2914 

-1.0856 

Swinoujscie(Pol) 

0.31731 

10.584 

6.5398 

5.8139 

Ustka(Pol) 

0.22654 

-0.01475 

0.080303 

0.097364 

Vaasa(Fin) 

0.21154 

0.118 

0.15485 

0.16146 

Visby(Swe) 

0.42269 

-0.00975 

0.16061 

0.19118 

W amemtinde(Ger) 

71.218 

-24.676 

13.101 

19.881 
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Table 6-5: Results of w 0 computation based on the orthometric heights of case 5. The results from one sta¬ 
tion to anther are very much different and therefore not acceptable. 


Station 

Name 

w 0 

(m 2 s~ 2 )/y 

1993.4-1990.8 

w 0 

(m 2 s- 2 )/y 

1997.4-1993.4 

w 0 

{m 2 s- 2 )/y 

1997.4-1990.8 

*0 

(mV 2 )/y 

mean 

Borkum(Ger) 

-1.5477 

-0.06875 

-0.65136 

-0.75594 

Degerby(Fin) 

0.12808 

2.8695 

1.7895 

1.5957 

Furuogrund(Swe) 

0.29846 

0.0835 

0.16818 

0.18338 

Hamina(Fin) 

0.075385 

-0.03175 

0.010455 

0.01803 

Hanko(Fin) 

0.10962 

8.2208 

5.0255 

4.4519 

Helgoland(Ger) 

-0.0038462 

0.11775 

0.069848 

0.061251 

Helsinki(Fin) 

0.19615 

0.94775 

0.65167 

0.59852 

Kemi(Fin) 

0.20423 

14.317 

8.7574 

7.7596 

Klagshamn(Swe) 

0.24923 

0.0955 

0.15606 

0.16693 

Klaipeda(Lit) 

- 

- 

- 

- 

Kronstadt(Rus) 

- 

- 

- 

- 

List/Sylt(Ger) 

-0.24923 

-0.03425 

-0.11894 

-0.13414 

Mantyluoto(Fin) 

- 

- 

- 

- 

Molas(Lit) 

- 

0.03675 

- 

- 

OlandsN.U.(Swe) 

0.21115 

0.027 

0.099545 

0.11257 

Raahe(Fin) 

-1.0954 

3.6047 

1.7532 

1.4208 

Ratan(Swe) 

-0.36615 

0.076 

-0.098182 

-0.12945 

Shepelevo(Rus) 

- 

- 

- 

- 

Spikama(Swe) 

0.17 

0.0615 

0.10424 

0.11191 

Stockholm(Swe) 

0.47231 

-2.4375 

-1.2912 

-1.0855 

Swinoujscie(Pol) 

0.31731 

10.584 

6.5398 

5.8139 

Ustka(Pol) 

0.22654 

-0.01475 

0.080303 

0.097364 

Vaasa(Fin) 

0.21154 

0.11775 

0.1547 

0.16133 

Visby(Swe) 

0.42269 

-0.00975 

0.16061 

0.19118 

W amemiinde(Ger) 

71.218 

-24.676 

13.101 

19.881 


7. Height Datum Difference of the 
Countries Around the Baltic Sea 

Indeed, the results we have presented in Table 5-1 
provide us with a means to estimate the difference 
between the height datums of countries around the 
Baltic Sea. If now we rearrange Table 5-1 in a way 
that the w 0 values of the each country be kept to¬ 
gether, we arrive at Table 7-1 , which helps us to 
see the variation of w 0 values from one country to 
the other. As it is shown in Table 7-1 for each 
country we have computed a reference potential by 
averaging the two closest w 0 at tide gauge stations 
of that country. Based on those reference potential 


values we computed the height datum difference 
between Finland, Germany, Lithuania and Sweden 
as shown in Table 7-2 in gravity space in terms of 
potential difference, and in Table 7-3 in geometry 
space, in terms of height difference. To convert the 
potential differences into height differences (i.e. to 
transformation from gravity space to geometry 
space) we use the mean value of vertical gradient of 
gravity -9.81802523 m / s 2 in Baltic Sea area, 

based on ellipsoidal harmonic model of de¬ 
gree/order 360/360. We leave the details of com¬ 
putation of vertical gradient of gravity to the next 
section. 



48 


A. Ardalan et al. 


Table 7-1: Gravitational potential U(\<t>,u), centrifugal potential V(X,</>,u), free-air gravity potential 
A W (A, </>, u]H 0 ), computed based on orthometric heights of the case 1 (Af. Poutanen et al ., 1999), and gauge 
value of geoid potential W 0 for the epoch of 1997.4. The highlighted stations are those used for computing 
the reference potentials. 


Station 

Name 

Gravitational 

Potential 

U(X^u) 

(m 2 s“ a ) 

Centrifugal 

Potential 

V F (A, it) 

(m 2 s~ 2 ) 

Free Air 

Reduction 

6W{\&u) 

(™ 2 s a ) 

Gauge 

Value 

Wt 

(mV*) 

Borkum (Ger) 

62598480.049 

38330.947 

44.888 

62636855.884 

Helgoland (Ger) 

62599593.036 

37218.924 

44.468 

62636856.427 

List/Sylt (Ger) 

62601102.613 

35713.291 

40.782 

62636856.686 

Wamemunde (Ger) 

62599530.263 

37213.313 

111.089 

62636854.665 



mean(m 2 s 2 ) 

62636856.557 




std(m 3 s 2 ) 

0.129 

Degerby (Fin) 

62609699.381 

27125.102 

27.738 

62636852.221 

Hamina (Fin) 

62610584.135 

26255.426 

16.015 

62636855.576 

Hanko (Fin) 

62609334.840 

27467.913 

50.254 

62636853.007 

Helsinki (Fin) 

62609868.216 

26924.835 

63.039 

62636856.090 

Kemi (Fin) 

62618330.677 

18455.223 

69.665 

62636855.565 

Mantyluoto (Fin) 

62612227.421 

24604.328 

24.638 

62636856.388 

Raahe (Fin) 

62616878.076 

19941.164 

34.653 

62636853.893 

Vaasa (Fin) 

62614577.940 

22265.875 

11.589 

62636855.403 



mean(?n 2 s 2 ) 

62636855.515 




std(m 2 s 2 ) 

0.056 

Furuogrund (Swe) 

62617205.722 

19542.016 

107.773 

62636855.511 

Klagshamn (Swe) 

62602014.413 

34818.792 

20.004 

62636853.208 

Glands N. U. (Swe) 

62605214.649 

31601.729 

40.514 

62636856.892 

Ratan (Swe) 

62615930.933 

20910.715 

15.077 

62636856.725 

Spikama (Swe) 

62613440.933 

23395.355 

18.895 

62636855.183 

Stockholm (Swe) 

62608441.995 

28295.629 

116.888 

62636854.512 

Visby (Swe) 

62605699.404 

31135.656 

19.556 

62636854.616 



mean (m 2 s 2 ) 

62636856.809 




std(m 2 s 2 ) 

0.083 

Klaipeda (Lit) 

62602168.314 

34409.333 

276.885 

62636854.532 

Molas (Lit) 

62602356.570 

34452.784 

44.926 

62636854.279 



mean(m 2 ,s 2 ) 

62636854.406 




std(m 2 s~ 2 } 

0.126 

Kronstadt (Rus) 

62609616.854 

27194.595 

- 

— 

Shepelevo (Rus) 

62609584.472 

27231.134 

- 

- 

Swinoujscie (Pol) 

62599126.7% 

37663.497 

— 

— 

Ustka (Pol) 

62600356,210 

36479.268 

14.997 

62636850.475 
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Table 7-2: Datum difference between the countries around Baltic Sea, in potential units. 



Finland 

Germany 

Lithuania 

Sweden 

Finland 

0 

(-1.042 ± 0.070) 

(1.109 ± 0.069) 

(-1.294 ± 0.050) 



(mV 2 ) 

(mV 2 ) 

(mV 2 ) 

Germany 

(1.042 ± 0.070) 

0 

(2.151 ± 0.090) 

(-0.252 ± 0.077) 


(mV 2 ) 


(mV 2 ) 

(mV 2 ) 

Lithuania 

(-1.109 ± 0.069) 

(-2.151 ± 0.090) 

0 

(-2.403 ± 0.075) 


(mV 2 ) 

( mV 2 ) 


(mV 2 ) 

Sweden 

(1.294 ± 0.050) 

(0.252 ± 0.077) 

(2.403 ± 0.075) 

0 


(mV 2 ) 

(mV 2 ) 

(mV 2 ) 



Table 7-3: Datum difference between the countries around Baltic Sea, in length unit, based on mean value of 
vertical gradient of gravity -9.81802523 m / s 2 in Baltic Sea area. 



Finland 

Germany 

Lithuania 

Sweden 

Finland 

0 

(-0.106 ± 0.007) 

(m) 

(0.113 ± 0.007) 

(m) 

(-0.132 ± 0.005) 

(m) 

Germany 

(0.106 ± 0.007) 
(m) 

0 

(0.219 ± 0.009) 

(m) 

(-0.026 ± 0.008) 

(m) 

Lithuania 

(-0.113 ± 0.007) 
(m) 

(-0.219 ± 0.009) 
(m) 

0 

(-0.245 ± 0.008) 

(m) 

Sweden 

(0.132 ± 0.005) 
(m) 

(0.026 ± 0.008) 
(m) 

(0.245 ± 0.008) 

(m) 

0 


8. Sea Surface Topography Map 
of Baltic Sea 

Having accessed to the MSL information at various 
tide gauge stations around Baltic Sea we have op¬ 
portunity to derive the Sea Surface Topography 
(SST) of the Baltic Sea. Here we use the GPS ob¬ 
servations of the Baltic Sea Level project, 3 rd cam¬ 
paign, which are the most accurate ones, and the 
orthometric heights of the case 1, which are the di¬ 
rectly observed heights above MSL in height datum 
of various countries (see Table 5-7). Of course now 
that we have estimated the datum difference be¬ 


tween the countries around Baltic Sea we can unify 
the datum of the orthometric heights of the case 1. 

The difference between the w 0 values presented 
in Table 5-1 and the average value w 0 = 
(62636855.75 ±0.21) (m 2 /s 2 ), provide us with 

the apparent SST, i.e., the deviation of the sea sur¬ 
face from the geoid in various tide gauge stations 
around the Baltic Sea, plus the height datum differ¬ 
ence. These potential deviations can be converted 
into metric units according to the transformation 
relation outlined in Table 8-1. 
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Table 8-1: Transformation relation of potential difference into height difference. 


“decomposition of actual geoid potential value w 0 
into the apparent geoid potential value at the tide gauges W 0i , 
and the disturbing part SW ” 


w 0 = + SW 

= w 0i + • k - u) 

__ _ 1 dW , , , 

= w 0i + • ^g^{u Q - u ) 


= Wn + 


du 

1 dW 


• Au® 


A u^ 


■\JQuu du 
the physical height difference 
(i.e. SSTT 
Oo - W 0i )J~g^; 


dW 

du 

“subject to” 


dW 

du 


■■ ^ Y2 Unm 


9q: h q 


Qn\m\i. g) 

+ufuy]u 2 4- £ 2 cos 2 (j) 


du 


,(A,0) 


( 8 . 1 ) 


( 8 . 2 ) 


(8.3) 


As shown in Table 8-1 the geoid potential value 
w 0 can be decomposed into the computed (or ap¬ 
parent) geoid potential value at the various tide 
gauges W 0i (index % runs from one to the total 
number of tide gauge stations), and the disturbing 
part SW . The disturbing part is due to the Sea 
Surface Topography (u () — u). The Sea Surface 
Topography (u 0 — u) is in terms of Jacobi ellipsoi¬ 
dal coordinates, and as such is not a physical com¬ 
ponent. However, A^ (1) = -sfg^(uQ — u) , derived 
by using the directional derivative operator 
(see(8.1) - (8.2)), is the required physical height. 
For the derivative of the potential with respect to 
u appearing in (8.2), we use the high-resolution 
model of (8.3). 

Table 8-2 presents the computed SST, A^ (1) , of 
the Baltic Sea at various tide gauge stations. We 


have also corrected the SST of tide gauge stations 
for the difference between the national height da¬ 
tums, to reach to the last column of Table 8-2, 
which shows the SST of the Baltic Sea Level tide 
gauges in the German height datum. Note that the 
SST could be more stable with respect to the Ger¬ 
man height datum than the height datums of the 
countries at the northern part of the Baltic Sea, for 
example, especially due to the runoff of rivers in 
the north in the spring period. That is why we have 
computed the SST of the Baltic Sea in the height 
datum of Germany. A contour map plot of com¬ 
puted SST of the Baltic Sea is shown in Figure 8-1. 

Table 8-3 shows the comparison of the com¬ 
puted SST for Baltic Sea as explained above and 
the one computed by J. Kakkuri M. Poutanen 
(1997). The difference between the two SST solu¬ 
tions is (0.011 d= 0.054) m . 
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Table 8-2: Vertical gradient of Potential, difference between apparent W Qi based on height of the stations in 


their respective national height systems and average value of geoid potential w 0 , apparent SST at the tide 
gauge stations of Baltic Sea level project in different height datums, corrections to German’s height datum, 
corrected SST in German height datum. 

Station 

Name 

1 dW 

- J~gm du 
(■ m / s 2 ) 

(wo-W 0 ) 

(m 2 /s 2 ) 

Aw' 1 * 

( w 0 -W 0 ) 

1 dW 

■yj CjuU dll 

(m) 

Correction to 
German height 
datum 
(m) 

Corrected 

SST 

(m) 

Borkum(Ger) 

-9.81363626 

-0.134 

0.0137 

0 

-0.014 

Degerby(Fin) 

-9.81862941 

3.529 

-0.3594 

0.106 

0.253 

Fumogrund(Swe) 

-9.82251837 

0.239 

-0.0243 

-0.026 

0.050 

Hamina(Fin) 

-9.81923277 

0.174 

-0.0177 

0.106 

-0.088 

Hanko(Fin) 

-9.81899359 

2.743 

-0.2794 

0.106 

0.173 

Helgoland(Ger) 

-9.81407564 

-0.677 

0.0690 

0 

-0.069 

Helsinki(Fin) 

-9.81916005 

-0.340 

0.0346 

0.106 

-0.141 

Kemi(Fin) 

-9.82306382 

0.185 

-0.0188 

0.106 

-0.087 

Klagshamn(Swe) 

-9.81537974 

2.542 

-0.2590 

-0.026 

0.285 

Klaipeda(Lit) 

-9.81547717 

1.218 

-0.1241 

0.219 

-0.095 

Kronstadt(Rus) 

-9.81913536 

- 

- 

- 

- 

List/Sylt(Ger) 

-9.81521803 

-0.936 

0.0954 

0 

-0.095 

Mantyluoto(Fin) 

-9.82002342 

-0.638 

0.0650 

0.106 

-0.171 

Molas(Lit) 

-9.81549953 

1.471 

-0.1499 

0.219 

-0.069 

OlandsN.U.(Swe) 

-9.81683093 

-1.142 

0.1163 

-0.026 

-0.090 

Raahe(Fin) 

-9.82228525 

1.857 

-0.1891 

0.106 

0.083 

Ratan(Swe) 

-9.82207194 

-0.975 

0.0993 

-0.026 

-0.073 

Shepelevo(Rus) 

-9.81897667 

- 



- 

Spikama(Swe) 

-9.82068805 

0.567 

0.0577 

0.026 

0.084 

Stockholm(Swe) 

-9.81840031 

1.238 

-0.1261 

-0.026 

0.152 

Swinoujscie(Pol) 

-9.81406739 

- 

- 

- 

- 

Ustka(Pol) 

-9.81472950 

5.275 

-0.5375 



Vaasa(Fin) 

-9.82096810 

0.347 

-0.0353 

0.070 

07)71 

Visby(Swe) 

-9.81722771 

1.134 

-0.1155 

-0.026 

0.142 

Wamemunde(Ger) 

-9.81434182 

1.085 

-0.1106 

0 

0.111 


Table 8-3: Comparison of our SST with the SST computed by J . Kakkuri M. Poutanen (1997) for the Baltic 
Sea. 


Station 

Name 

our SST 
(m) 

SST (/. Kakkuri 

M. Poutanen , 
1997) 

Difference 
between 
two SST’s 

Borkum(Ger) 

-0.014 

- 

- 

Degerby(Fin) 

0.253 

0.095 

0.158 

Furuogrund(Swe) 

0.050 

0.054 

-0.004 

F[amina(Fin) 

-0.088 

0.052 

-0.140 

Hanko(Fin) 

0.173 

0.079 

0.094 

Helgoland(Ger) 

-0.069 

- 

- 

Helsinki(Fin) 

-0.141 

0.024 

-0.165 

Kemi(Fin) 

-0.087 

0.146 

-0.233 

Klagshamn(Swe) 

0.285 

-0.363 

0.648 

Klaipeda(Lit) 

-0.095 

- 

- 
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Tabelle 8-3: (continued) 


Station 

Name 

our SST 

(m) 

SST (7. Kakkuri 

M. Poutanen , 
1997) 

Difference 
between 
two SST’s 

Kronstadt(Rus) 

- 

- 

- 

List/Sylt(Ger) 

-0.095 

- 

- 

Mantyluoto(Fin) 

-0.171 

0.044 

-0.215 

Molas(Lit) 

-0.069 

- 

- 

OlandsN.U.(Swe) 

-0.090 

-0.03 

-0.060 

Raahe(Fin) 

0.083 

0.067 

0.016 

Ratan(Swe) 

-0.073 

0.041 

-0.114 

Shepelevo(Rus) 

- 

0.260 

- 

Spikama(Swe) 

0.084 

0.116 

-0.032 

Stockholm(Swe) 

0.152 

0.199 

-0.047 

Swinoujscie(Pol) 

- 

-0.300 


Ustka(Pol) 

- 

-0.344 


Vaasa(Fin) 

-0.071 

-0.01 

0.061 

Visby(Swe) 

0.142 

0.050 

0.092 

W amemiinde(Ger) 

0.111 

-0.143 

0.254 



Figure 8-1: The computed SST map of Baltic Sea in German Height Datum. Equi¬ 
distant Conic Map Projection with two standard parallels: 55 ° N and 65 ° N , and 
reference ellipsoid of WGD2000. 
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9. Conclusions 

Let us summaries the results we obtained in various 
sections as follows. 

(i) Our best estimate for w 0 value based on the 
GPS observations of the Baltic Sea Level 
project, 1 st , 2 nd , and 3 rd campaigns, and tide 
gauge observations is 

w 0 = (62636855.75 ± 0.21 )(m 2 /s 2 ) 

(ii) Our best estimate for w 0 value based on the 
GPS observations of the Baltic Sea Level 
project, 1 st , 2 nd , and 3 rd campaigns, and tide 
gauge observations is 

Wq = (-0.0099 ± 0.00079) (ra 2 /s 2 )/year 
or 

w 0 / 7 = 1 . 0 (mm / year) 


(iii) Amongst the different geoid solutions 
proposed for Baltic Sea the one introduced 
in M. Poutanen et al. (1999) is the most 
accurate one. 

(iv) The geoid proposed by J. Kakkuri (2000) 
has a shift of approximately 0.244 m . 

(v) As the comparison with already published 
w 0 values, we refer Table 9-1. 

Finally, based on the results obtain we can conclude 
that our method is quite successful for the compu¬ 
tation of fundamental geodetic parameter w 0 , the 
potential value of Gauss-Listing geoid, as well as 
its time derivative w 0 . Besides, the proposed meth¬ 
odology can also be quite helpful in (i) unification 
of the National Height Datums (ii) computation of 
the high-resolution Sea Surface Topography maps , 
and (iii) accuracy estimation of the geoid solutions 
tailored to Sea areas. 


Table 9-1: Comparison of potential value of geoid w 0 computed by various authors. 


Author 

w 0 value ( m 2 /s 2 ) 

Data source / Computation method 

D. Nesvorny,and Z. 
Sima (1994) 

62636857.5 ± 1.0 

satellite altimetry data 

Bursa et al. (1997a) 

62636855.72 ± 0.5 

satellite altimetry data, gauge stations values 

Bursa et al. (1997b) 

62 636 855.80 ± 0.5 

satellite altimetry data, gauge stations values 

Grafarend and Ar¬ 
dalan (1997) 

62636855.8 zb 3.6 

ellipsoidal harmonic expansion and, and tide gauge 
information of Baltic Sea, GPS observations of BSL 
projects 2 nd campaign 

Bursa et al. (1998) 

62 636 855.611 ± 0.5 

TOPEX/POSEIDON altimeter data 

Bursa et al. (2000) 

62636856.0 ± 0.5 

TOPEX/POSEIDON altimeter data 

Our results 

62636855.75 ± 0.21 

ellipsoidal harmonic expansion, and tide gauge in¬ 
formation of Baltic Sea, GPS observations of BSL 
projects 1 st 2 nd and 3 rd campaigns 
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Appendix A: 

Spheroidal Coordinates 

It has been revealed by the Great Early 18 th Cen¬ 
tury Expeditions that the Earth is not geometri¬ 
cally a sphere, but nearly an oblate ellipsoid-of- 
revolution >a b . (See for example J. Kakkuri et 

al (1986), J.R. Smith (1986) and E. Tobe (1986) 
as the historical review of the progress in deter¬ 
mination of the shape of the Earth). Therefore, 
representation of gravity field of the Earth in 
terms of ellipsoidal harmonics is more accurate 
and even more convenient than for example 
spherical coordinates and spherical harmonics. 

Here briefly we will review the main features of 
ellipsoidal coordinates and ellipsoidal harmonic 
expansion and invite the interested readers for 
more details to a visit to N. Thong and W. Gra¬ 
farend (1989). 


Definition A-l: (spheroidal coordinates {A, 4>,u}) 


In terms of ellipsoidal coordinates {A, </>, u) , a point in space can be located as the intersection point of the 
following family of surfaces. 


(i) the family of confocal, oblate spheroids 


E : 


V« 2 +£ 2 ,u 


:= xel 


u 2 +£ 2 


H —% = 1, u € ( 0 ,+ oo), e 2 := a 2 — b 2 


(A.l) 


(ii) the family of confocal half hyperboloids 


cos 0,6 sin 0 


= xel 


x 2 + y 2 

£ 2 COS 2 (p 


Z , , , 7T 7T, 

- X<t> e h 2 ’ 2 1 ’ * 


(Hi) the family of half planes 

!Pc 2 osA,smA := { xeM 3 I y = a; tan A, A e [0, 2 tt]} 


(A.2) 


(A.3) 


As shown in Figure A-l longitude A gives orien¬ 
tation to the half planes, latitude f is the inclina¬ 
tion of asymptotes of the confocal half hyperbo¬ 


loids, and the elliptic coordinate u is the semi-minor 
axis of confocal oblate spheroids (confocal, oblate 
ellipsoids of revolution). 
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Table A-l: Forward and backward transformation ofCartesian coordinates {x,y,z} into ellipsoidal coordi¬ 
nates {\cj),u} 


(i) Forward transformation 

{\<I>M {x,y,z} 

x = yju 2 -he 2 cos(/>cosA 

y = \lu 2 + e 2 cos 0 sin A 
z = usincj) 


(ii) Backward transformation of { x , y, z } {A, </>, u} 

y 

arctan—for x > 0 and y > 0 
x 

V 

arctan — + n for x < 0 and y ^ 0 
x 

V 

A = <1 arctan — + 2tt for x > 0 and y < 0 
x 

— for x = 0 and y > 0 

2 y 

3 ^ for x = 0 and y < 0 


4> = (sgn z) arcsin | ^ 


[e 2 - (x 2 + 2/ 2 + * 2 ) 
1/2 


+ V (^ 2 + </ 2 +* 2 -£ 2 ) 2 + 4eV]} 


ti = e * 2 +y 2 +z 2 -s 2 


-i 


+sl{x 2 +y 2 +z 2 


t 2 f +4 eV]} 1/2 


(A.4) 


(A.5) 


(A. 6 ) 

(A.7) 


Definition A-2: Basic geometry of ellipsoidal coordinates {A, o. //} 

(7) Jacobi matrix of the forward transformation {\,<f>,u} >—► {a.% j/. 2 } 

From equation (A.4) Jacobi matrix “ J ” of the transformation from ellipsoidal coordinates {A, o. « } into 
Cartesian coordinates { 2 ;, y, z] can be obtained as follows. 


J := 


A, 

Yx 

Z x 


X* 


Xu 

Y u 

Yu 


The partial derivatives used in (A.8) are defined as 

X x = D x X = —V u 2 + e 2 cos </> sin A 

Y x = D X Y = Vm 2 + e 2 cos <p cos A 
= 0 

= D^X = —\lu 2 + £ 2 sin <f> cos A 

= — Vu 2 + e 2 sin f sin A 
= D^Z = ucosf 


(A. 8 ) 



National Height Datum, the Gauss-Listing Geoid Level Value 57 


X, - D U X = 
Y u = D U Y = 


V u 2 + t 

U 


r cos p cos A 


cos p sin A 


G := J J 


( U 2 + £ 2 )cOS 2 p 

0 

0 


Ju 2 + £ 2 

= sin</>. 

(ii) The metric tensor 

\d\\ 

dS 2 = [d\,d(j),du\J*J\d(j) 
du 

0 0 

u 2 + e 2 sin 2 p 0 

0 ( u 2 + e 2 sin 2 0) / (w 2 + £ 2 ) 


(A.9) 


:= g n m Vn,ra = 1,2,3 (A. 10) 


(iii) Laplacian 


A = ±l± ( 2lljL + jL(giJL + JL(JlJL) 

•Jg 1 d\ ' g n d\ d<t) { g 22 d<t> du { g 33 du’ 

1 

--!-:-t.fl.n tt, 

dp 


u 2 + £ 2 sin 2 p[(u 2 + e 2 ) cos 2 p d\ 2 

d 2 0 d , 2 2 , a 2 ' 

+ 777T + 2tZ —-h (u + £ )-^~2 


(A.11) 


dp 2 


du 


du 2 


z 



Figure A-l: Spheroidal Coordinates {A ,p,u}. 

The coordinate surfaces are of the type: (i) sphe¬ 
roids (u = const.), (ii) half hyperboloids of one 
sheet (p = cost.) and, (iii) half planes (A = 

const.) 
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Appendix B: 

Normalised Associated Legendre 
Functions of First and Second Kind 


Here we define the normalised associated Legendre 
functions of the first kind P n * m (sin 0) by means of 
recurrence relations as follows. 


Pnn (sin (f>) = 1 cos <j> Cl (sin cj>) 


^ n *n-i(sm0) = ^ 2n + -* cos <j> (sin <f>) 


P^m (sin <j>) — 


V 4 n 2 — 1 


J n 2 - m. 


sin 4> Pn-i,m (sin <t>) 


<7(2 n + l)(n + m — l)(ra — m — 1) 
A /(n 2 "^-"m 2 7(2n"^-"3) 

subject to 

Vn G [3,oo) and m € [0,n — 2] 


Pn-2,m (sin <; 


(B.l) 

(B.2) 


(B.3) 


with starting values 

Poo (sin <f>) = 1 
Pi'o (sincj)) = V3sin </> 

Pi X ( sin</>) == V3cos</> 

P 2 o(sin<p) = ^(3sin 2 <j> - 1) 
P 2 i (sin <p) = Vl5 sin (p cos </> 

., . . Vl5 2 , 

F 22 (sin0) = ~ cos <j> 


(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.9) 


The associated Legendre functions of the second kind can be defined by an integral relation of the type 

(B.10) 


Ql(|) = i n+1 Q„m(i |) 


Qn 


u = (— l) m 2 n (n + m)\m \ , u 2 +£ 2 , m/2 _ °r 
e i n+1 (n — m)!( 2 m)! e 2 J 


sinh 2m r dr 


yju 2 + £ 2 


(B.l 1) 


0 (- + 

£ £ 


coshr 


\n+m+1 


However, in practice instead of the above integral formulas the associated Legendre functions of the second 
kind are calculated via the recursive relations which enjoy the numerical stable, especially for the higher de¬ 
grees and orders (N. Thong and E. Grafarend , 1989, G. Sona, 1996) 


Qn m© = 

6 k=0 


(B. 12) 


o* (u\ = g2(1 ~ ra ~ I m 1 ~ 2k )( n + I m 1 + 2fc ) Q* (u) w fc > 1 

Qn\ m \ k ( ) 2k(2n + 2k + l)u 2 > 


(B. 13) 
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Q*n\m\o( u ) = cosh H 7 f~) n+1 VnG#, me[-n,n] C Z (B.14) 

1 1 u 

The summation (B.12) is continued until 

^n|m|A:max (^) ^ n \ m \ k max-1 ^ ^ (B.15) 

with starting values for n- 0, 1,2 and m = 0 

Qo{~) = arccot(-) (B.16) 

£ £ 

Ql (-) = 1 - - arc cot (-) (B. 17) 

£ £ £ 

&*(-) - + l)arccot(-) - 3-] (B.18) 

£ 2 £ £ £ 


In (B.15) a can be selected according to required accuracy. In our calculations double precision accuracy, 

i.e., cr =1E-16 was adopted. 


Appendix C: Transformation 
of Spherical Harmonic Coefficients 
into Spheroidal/Ellipsoidal 
Harmonic Coefficients 

Nowadays it is a common practice to represent the 
“Standard Gravity Earth Models ” in terms of 
spherical harmonics. Fortunately, precise transfor¬ 
mation relations between spherical and ellipsoidal 
harmonic coefficients are available and therefore 
one can transfer the spherical harmonic coefficients 
into ellipsoidal ones without any loss of accuracy. 
Lemma C-l offers a summary of the transformation 
of spherical harmonic coefficients into ellipsoidal 
harmonic coefficients according to C. Jekeli (1981, 
1988). In conjunction with the ellipsoidal harmon¬ 
ics, contribution made by D. Gleason (1988, 1989), 
G. Sona (1996) and J. Yu and H. Cao (1996) can 
also be acknowledged. 

Lemma C-l: (Transformation of spherical har¬ 
monic coefficients into ellipsoidal harmonic coeffi¬ 
cients) 


Spherical harmonic coefficients, u n , ™ ( spherial ), 
can be uniquely transformed into ellipsoidal har¬ 
monic coefficients, u n , m (ellipsoid) Ql )H (-) are 

associated Legendre functions of the second kind , 
see equation (B.12) for the recursive relation. 

As one can read from (C.1)-(C.5) any ellipsoidal 
harmonic coefficient, u n , m {spheroid) is equal to 

the spherical harmonic coefficient of the same de¬ 
gree and order u n , m {sphere) plus a linear combi¬ 
nation of spherical harmonic coefficients of the 
lower degree but the same order. 

For details on transformation of spherical har¬ 
monic coefficients into ellipsoidal ones, we refer to 
A. Ardalan and E. Grafarend (2000). The ellipsoi¬ 
dal harmonic coefficients in different permanent 
tide systems can be accessed also from 
http://www.uni-stuttgart.de/gi/research/index.html# 
Projects. 


£ {n—m) / 2 

u n , m{ellipsoidal ) = Q* w (-) ^ \ K , w„_ 2W> ^{spherical) 

£ 1=0 


(2n — 21 )! n ! 

(2 n — 4/ + l)(n — m){n + m)! 

(2 n) ! l\ {n — 1)! 

(2n -K 1 ){n — 21 + m)\{n — 2 1 — m)\ 


Vn G [0,oo) and m G [— n, + n] 


(C.l) 

(C.2) 

(C.3) 
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By expanding the factorials in (C.2), one can reach to following recursive formula, which is numerically sta¬ 
ble especially for high degrees and orders n/m. 




((2 n — 4/ + 1 ){n — 21 — m 4- l)(n — 2/ — m + 2 )(n — 21 4- m + 1 )(n — 21 +m + 2)) 

2k(2n - 21 + 1)(2n - 4/ + 5) 1/2 
V l G [1, (n — m) / 2], n G [0, oo) and m G [—n, + n] 


1/2 


C £ K 


(C.4) 


with the start value 

\, m ,o =1 V n, m 


(C.5) 



